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ABSTRACT 


A  combination  of  creeping -wave  analysis  and  diffraction  theory 
has  been  developed  for  determining  the  radar  cross  section  of  bodies 
for  which  exact  solutions  are  not  available.  The  known  solutions  for 
the  perfectly  conducting  cylinder  and  sphere  have  been  used  to  specify 
attenuation  and  diffraction  coefficients  for  the  creeping  wave.  The 
creeping  wave  contribution  is  added  to  the  geometrical  optics  or 
physical  optics  contribution  from  the  specular  point  to  determine  the 
total  scattered  field.  It  is  demonstrated  that  this  type  of  solution 
is  applicable  to  ogives,  ogives  with  spherical  caps,  and  prolate 
spheroids.  Wedge  diffraction  theory  has  been  combined  with  creeping 
wave  analysis  to  calculate  the  edge-onbackscatter  of  circular  and 
ogival  disks.  It  is  necessary  to  modify  the  magnitudes,  but  not  the 
forms,  of  the  creeping  wave  attenuation  and  diffraction  coefficients 
when  treating  a  creeping  wave  on  an  edge.  This  form  of  analysis 
is  adaptable  to  calculation  of  the  radar  cross  section  of  composite 
bodies  where  both  volumetric  shapes  and  edges  may  occur  in  com¬ 
bination.  Once  the  specular  points,  wedge  diffraction  points,  and 
the  attachement  points  and  paths  of  the  creeping  waves  are  determined 
the  computation  of  the  scattered  field  is  straightforward. 
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MEMORANDUM  ON  ANALYSIS  OF  ECHO  AREA  OF  TARGETS 
USING  GEOMETRICAL  THEORY  OF  DIFFRACTION 
AND  CREEPING  WAVE  THEORY 

INTRODUCTION 

This  report  is  a  preliminary  attempt  to  extend  the  creeping-wave 
analysis  to  bodies  for  which  it  has  not  previously  been  applied.  Several 
new  approaches  incorporating  geometric  diffraction  theory  and  wedge 
diffraction  concepts  are  proposed  and  evaluated.  This  report  is  intended 
to  summarize  the  work  that  has  been  completed  to  date  and  to  indicate 
problem  areas  which  must  be  investigated  in  order  to  develop  a  consist¬ 
ent  c reeping-wave  theory  of  scattering.  It  is  not  intended  to  represent 
a  complete  analysis  of  the  problem. 

The  determination  of  the  radar  cross  section  of  targets  composed 
of  volumetric  shapes,  ducts,  fins,  and  edges  is  a  difficult  but  important 
problem.  One  approach  to  this  problem  is  to  approximate  the  target  by 
a  collection  of  shapes  for  which  the  individual  scattering  properties  are 
known  and  thus  estimate  the  scattered  field  in  a  given  direction  by 
locating  the  dominant  scatterers  and  adding  their  contributions. 

A  major  flaw  in  this  approach  is  that  a  general  technique  for 
analyzing  the  part  of  the  body  that  is  in  the  shadow  region  or  the  region 
not  directly  illuminated  by  the  incident  wave  is  lacking.  Peters1  showed 
that  a  traveling-wave  technique  is  applicable  for  long  thin  bodies  such  as 
the  ogive  where  the  dominant  scattering  mechanism  is  chiefly  influenced 
by  the  part  of  the  body  in  the  shadow  region.  This  technique  has  been 
used  by  all  researchers2  for  this  type  of  target  despite  objections  that 
it  lacks  mathematical  rigor.  Franz  and  Depperman,  Keller  and  Levy  * 
and  others  have  developed  the  creeping  wave  for  successfully  treating 
the  shadow  region  of  the  cylinder.  Similar  treatment  of  the  sphere  has 
recently  been  developed.6  However  these  solutions  are  restricted  to 
bodies  for  which  canonical  solutions  are  available.  The  major  purpose 
of  this  report  is  to  show  that  the  shadow  region  of  a  target  can  be  treated 
by  a  relatively  simple  procedure. 

Another  approach  is  the  use  of  point-matching  boundary  value 
techniques  to  generate  a  system  of  linear  equations  which  may  be  solved 
for  the  surface  currents  using  a  computer.  The  scattered  field  is  then 
computed  from  the  surface  currents.  In  this  second  approach  one  is 
limited  by  the  capacity  of  the  computer  and  the  long  computation  time. 
This  report  demonstrates  that  the  techniques  of  the  geometrical  theory 
of  diffraction  and  creeping  wave  theory  offer  a  means  by  which  the  scat¬ 
tering  of  a  composite  target  can  be  evaluated  to  engineering  accuracy 
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using  relatively  simple  formulations  which  do  not  require  large  amounts 
of  computation  time.  In  addition,  this  analysis  allows  the  evaluation  of 
interactions  between  the  component  parts  of  the  scatterer.  At  the  pre¬ 
sent  time  only  several  simple  scatterers  have  been  treated.  These  in¬ 
clude  the  ogive,  ,7  the  spherically  capped  ogive,  8  the  finite  cone,  10  the 

Q 

cone-sphere,  and  the  disk  for  edge-on  aspects. 

The  techniques  of  the  geometrical  theory  of  diffraction  and  creep¬ 
ing-wave  analysis  depend  upon  the  assumption  that  the  propagation  paths 
of  the  diffracted  and  creeping-wave  fields  may  be  described  using  ray 
tracing  and  that  diffraction  is  a  local  phenomenon.  These  assumptions 
make  it  possible  to  determine  the  contributing  scattering  mechanisms 
of  the  target  for  a  given  direction  and  also  allow  the  calculation  of  inter¬ 
actions  between  parts  of  the  scatterer.  The  scattered  fields  are  then 
obtained  by  summation  of  the  individual  contributions  . 

In  order  to  develop  useful  diffraction  and  attenuation  coefficients 
for  bodies  of  revolution,  the  scattering  of  the  perfectly  conducting  sphere 
has  been  evaluated  in  the  form  of  a  specular  point  contribution  plus  a 
c reeping- wave  contribution.  Physical  optics  has  been  assumed  to  apply 
to  the  specular  point  and  the  creeping-wave  diffraction  and  attenuation 
coefficients  have  been  obtained  by  matching  the  exact  back  scattered 
field  at  a  single  value  of  kQa.  These  values  were  then  used  to  calculate 
both  backscattered  and  bistatic  fields  for  a  range  of  kQa  values.  The 
excellent  result  obtained  indicated  that  such  a  format  was  useful.  This 
form  of  solution  was  then  extended  to  calculate  the  back  scattered  fields 
of  a  prolate  spheroid  with  good  results.  The  application  of  this  form  of 
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solution  to  the  ogive  and  spherically  capped  ogive  had  been  done  previously  ’ 
and  is  in  agreement  with  the  formulation  proposed  here. 

The  impulse  technique  has  also  been  applied  to  obtain  the  scattered 
fields  for  the  prolate  spheroid^  and  for  the  cone  sphere.9  Basically  the 
impulse  solution  for  the  cone  sphere  is  forced  to  contain  the  creeping- 
wave  contribution  by  transforming  to  the  time  domain  the  exact  scattered 
fields  for  the  sphere,  after  subtraction  of  the  specular  scattered  fields. 

The  impulse  solution  for  the  prolate  spheroid  makes  use  of  the  required 
behavior  in  the  time  domain  of  the  scattered  fields  to  treat  the  shadow 
region  of  the  spheroid;  i.e.,  the  moment  conditions  are  enforced.  This 
is  a  clever  and  rewarding  technique  which  should  be  pursued  further. 

In  the  present  report  the  earlier  studies  applying  the  geometrical 
theory  of  diffraction  technique  to  these  non- canonical  shapes  are  sum¬ 
marized  and  certain  significent  points  are  clarified.  In  addition  certain 
new  mechanisms  are  introduced. 
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The  wedge  diffraction  coefficient  has  been  applied  to  the  problem 
of  determining  the  creeping-wave  fields  which  are  transmitted  and 
reflected  at  a  discontinuity  of  the  first  derivative  of  the  surface  on  a  body 
of  revolution.  Results  obtained  using  this  formulation  indicate  that,  to  a 
first  approximation,  the  wedge  diffraction  coefficient  is  applicable  to 
bodies  of  revolution,  and  thus  to  three-dimensional  scattering  in  spite 
of  the  two-dimensional  nature  of  the  wedge  scattering  problem.  The 
wedge  diffraction  coefficient  has  also  been  applied  in  combination  with 
creeping  waves  to  determine  the  scattered  fields  of  a  circular  disk  and 
a  "flat  ogive"  with  excellent  results.  These  results  indicate  that  an 
extension  of  these  techniques  to  the  analysis  of  finite  edges,  such  as 
occur  on  fins,  is  practical.  Furthermore,  the  relatively  simple  form 
of  the  equations  involved  in  these  computations  indicates  that  this  approach 
could  be  useful  in  identifying  the  geometry  in  the  shadow  region. 

The  creeping  wave  and  diffraction  form  of  analysis  permits  a  simple 
physical  picture  of  the  scattering  process.  The  calculation  of  the  scat¬ 
tered  field  is  straightforward,  with  the  diffraction  and  attenuation  coef 
ficients  being  simple  formulas  to  evaluate,  or  program  for  automatic 
computation.  In  general  the  time  required  to  compute  a  full  360°  back 
scatter  pattern  in  increments  of  2°  is  less  than  one  minute  on  an  IBM 
7094.  Thus  this  technique  seems  valuable  in  the  determination  of  radar 
cross  section  of  complex  targets* 

New  solutions  presented  in  this  report  include  the  flat  ogive, 
circular  disk  for  edge-on  incidence,  the  ogive  of  revolution,  and  the 
prolate  spheroid. 


I.  CREEPING  WAVE  ANALYSIS  OF 
BODIES  OF  REVOLUTION 

A .  Review  of  the  Analysis  of  a  Sphere 


The  concept  of  creeping  waves  was  introduced  by  Franz  and 
Depperman  9  for  the  interpretation  of  the  scalar  solution  for  diffraction 
by  a  circular  cylinder  or  sphere.  Senior  and  Goodrich14  have  obtained  a 
representation  similar  to  that  of  Franz  and  Depperman  through  the  appli¬ 
cation  of  the  Watson  transformation  to  the  Mie  series  solution  for  the 
sphere.  Moreland,  Peters,  and  Kilcoyne  have  applied  approximate  dif¬ 
fraction  and  attenuation  coefficients  developed  by  Peters11  to  calculate  both 
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monostatic  and  bistatic  cross  sections  of  the  sphere.  Kouyoumjian  has 
also  presented  a  creeping-wave  solution  for  the  sphere  which  includes  all 
higher-order  modes.  The  approach  presented  here  is  that  of  Moreland, 
Peters,  and  Kilcoyne. 
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Referring  to  Fig.  1  we  see  that,  according  to  the  geometrical  theory 
of  diffraction,  a  portion  of  the  energy  incident  upon  the  sphere  is  trapped 
at  the  point  of  tangency,  point  A.  The  ratio  of  trapped  field  to  incident 
field  is  defined  as  the  diffraction  coefficient.  The  trapped  energy  then 
travels  about  the  curved  body  and  is  attenuated  by  radiation  tangential 
to  the  body  as  it  travels.  Finally,  energy  is  radiated  along  a  tangent 
to  the  body  in  the  scattering  direction  of  interest.  The  diffraction  and 
attenuation  coefficients  which  describe  this  process  are  a  function  of 
the  radii  of  curvature  of  the  body  in  wavelengths.  The  diffracted  field  for 
the  lowest-order  creeping-wave  mode  is  given  by 


,1) 


Ud 

Ui 


-j 


DiDt  e 


fB 

k0L  +  ^  a{p)ds 
"A 


F(ri,  r2  )  , 


where 

Ui  is  the  incident  field, 

Ud  is  the  diffracted  field, 

Dj  is  the  diffraction  coefficient  at  the  incidence  point, 

Dt  is  the  diffraction  coefficient  at  the  point  of  reradiation, 
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&(p )  is  the  complex  attenuation  coefficient  as  a  function  of 
local  radius  of  curvature, 
ds  is  the  differential  arc  length, 

L  is  the  path  length  along  the  body, 
kQ  is  the  free  space  propagation  factor,  and 
F(ri,  r2)  is  the  phase  function  relating  source  and  observation 
points  • 

If  the  body  has  constant  curvature,  as  for  the  sphere,  the  radius  of 
curvature  at  points  A  and  B  are  equal  and  the  complex  attenuation  co¬ 
efficient  is  constant#  Thus  Eq.  (1)  can  be  written  as 


(2) 


Ud 

Ui 


=  D2  e"^-  jk°L  +  a  LJ 


F(ri,  r2)  . 


Cases  for  which  Eq .  (1)  must  be  used  are  illustrated  later  as  for  the 
prolate  and  oblate  spheroids.  In  such  a  case  the  local  radii  of  curvature 
at  points  A  and  B  must  be  used  to  evaluate  Di  and  Dt>  and  the  integral  in 
Eq.  (1)  must  be  evaluated. 

The  diffraction  and  attenuation  coefficients  have  been  derived  by 
comparing  solutions  in  this  form  to  canonical  solutions  for  cylinders, 
elliptic  cylinders,  and  spheres.5  This  simple  description  is  compli¬ 
cated,  however,  by  the  presence  of  higher-order  modes  and  the  existence 
of  higher-order  terms  in  D  and  Ot  which  depend  upon  the  orthogonal  radius 
of  curvature  . 

In  the  case  of  the  sphere,  Keller16  introduces  an  additional  factor 
to  account  for  the  convergence  of  rays  as  they  propagate  along  a  geodesic 
(great  circles  for  the  sphere).  Keller  has  evaluated  this  spreading  factor 
by  comparison  with  the  asymptotic  form  of  the  exact  solution  of  the  scat¬ 
tered  fields  by  the  sphere.  He  finds  that  the  diffraction  coefficient  and 
the  attenuation  coefficient  are  identical  to  those  for  the  cylinder.  He  also 
introduces  a  phase  jump  which  results  from  the  caustic  formed  at  the  rear 
of  the  sphere  by  the  rays  propagating  along  the  geodesic.  Senior  has 
extended  the  formulation  of  the  creeping-wave  component  of  the  fields 
scattered  by  a  sphere  of  small  radius  by  introducing  the  second  term  in 
the  asymptotic  expansion.  His  results  are  in  excellent  agreement  with 
the  exact  values,  but  his  analysis  cannot  be  written  in  a  form  utilizing 
diffraction  coefficients  . 

The  approach  suggested  by  Moreland,  Peters  and  Kilcoyne,8  and 
adopted  here,  differs  from  the  above,  in  that  the  exact  magnetic  fields  at 
the  surface  of  the  sphere  are  used  to  determine  the  complex  attenuation 
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coefficient.  This  approach  should  include  all  modes  if  an  accurate  approxi¬ 
mation  to  the  exact  magnetic  fields  on  the  surface  can  be  generated.  It 
has  been  found  that  throughout  the  resonance  region  the  attenuation  factor 
for  the  creeping  wave  is  approximately 


where 


a  is  the  sphere  radius, 

kQ  is  the  free  space  propagation  factor,  and 
QoC  is  the  attenuation  factor  for  the  lowest  order 
creeping-wave  mode  on  the  cylinder. 

The  diffraction  coefficient,  D,  has  been  found  by  multiplying  the 
diffraction  coefficient  of  the  lowest-order  creeping  wave  on  the  cylinder 
by  a  constant.  The  constant  is  determined  by  matching  the  exact  creep¬ 
ing-wave  fields  on  the  sphere  using  the  attenuation  coefficient  given  by 
Eq .  (3).  The  diffraction  coefficient  is  then  found  to  have  the  form 

TT 

2  J_  ^  1  / 

(4)  Ds  =  0.27  a3  X3  e 


The  process  through  which  these  values  for  the  diffraction  and 
attenuation  coefficients  were  obtained  is  simple  and  is  as  follows:  The 
magnetic  field  shown  in  Fig.  2  in  the  shadow  region  on  the  surface  of  the 
sphere  is  written  as  the  sum  of  two  creeping  waves  as 


(5) 


H(j)  =  H0 


(jk0  +  as)  -  ej 


+  e 


(iko  +  as)  a(j  +  e) 


and  the  radial  electric  field  may  be  written  as 
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Now  H0  and  E0  are  chosen  to  match  the  exact  solution  as  closely  as  possible 
for  some  value  of  kQa.  A  value  of  kpa  =  10  was  chosen.  Then,  assuming 
a  s  to  be  of  the  form  of  aoc  one  obtains  Eq.  (3).  In  slightly  abbreviated 

fashion  one  may  write 


.  T T 

i  i  Jr 

(7)  as  =  0.84  a'3  X  ‘3  e  6 


Next,  the  exact  scattered  fields  are  compared  to  the  fields  predicted  using 
Eq .  (2)  and  the  magnitude  of  D2  is  determined. 

The  form  of  Eqs  .  (5)  and  (6)  assumes  that  the  creeping  waves  are 
launched  at  points  0  =  it/  Z  and  0=  -tt/  2  and  travel  in  opposite  directions 
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about  the  sphere.  For  other  attachment  points,  the  paths  are  parallel  to 
this  great  circle  at  the  point  of  attachment.  The  direction  of  the  creeping 
waves  at  the  point  of  attachment  is  fixed  by  the  Poynting  vector  of  the 
incident  field.  After  the  wave  is  attached  to  the  surface  the  minor  creep¬ 
ing  wave  attenuates  rapidly,  while  the  major  creeping  wave  is  attenuated 
slowly.  Thus  the  resulting  surface  wave  at  some  point  removed  from  the 
shadow  boundary  has  only  a  radial  component  of  electric  field.  However 
the  existence  of  both  9  and  0  components  of  the  magnetic  field  at  the  point 
of  attachment  causes  a  tilt  of  the  Poynting  vector  on  the  surface.  For 
points  on  the  shadow  boundary  removed  from  the  major  wave  axis  this 
tilt  precludes  the  convergence  of  rays  to  form  a  single  caustic  at  the 
rear  of  the  sphere. 

This  approximate  form  may  be  compared  to  the  exact  solution  for 
the  fields  on  the  surface  of  the  sphere  along  this  great  circle  route 
through  straightforward  computation.  These  fields  are  illustrated  in 
Fig.  3  for  a  range  of  kQa  values.  It  is  seen  that  fields  of  the  form  given 
in  Eqs  .  (5)  and  (6)  are  a  reasonable  approximation  to  the  exact  fields. 

In  addition  the  backscattered  fields  due  to  the  creeping  wave  may  be 
written  as 


(8) 


where  R  is  the  observation  distance  from  the  center  of  the  sphere.  The 
magnitude  and  phase  of  the  creeping  wave  obtained  using  Eq .  (8)  can  be 
compared  to  the  magnitude  and  phase  of  the  exact  creeping-wave  fields 
obtained  by  subtracting  the  specular  point  contribution  of  the  physical 
optics  term  from  the  exact  scattered  fields;  i.e., 


(9) 


suppressed . 


The  magnitudes  of  the  exact  and  approximate  c reeping- wave  components 
are  compared  in  Fig.  4  and  their  phases  are  compared  in  Fig.  5.  The 
agreement  is  good  for  a  wide  range  of  kQa  values,  thus  establishing  the 
usefulness  of  Eq.  (8). 

The  concept  of  two  creeping  waves  propagating  around  the  sphere 
may  also  be  used  to  determine  the  bistatic  scattering  in  the  plane  of 
incidence  in  a  simple  manner.  Referring  to  Fig.  6  we  see  that  the  com¬ 
ponents  of  the  bistatic  field  are 
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Fig.  3.  Normalized  fields  at  the  surface  of  a  conducting  sphere. 
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Fig.  3.  Normalized  fields  at  the  surface  of  a  conducting  sphere. 
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ESp,  the  specularly  reflected  field  (where  it  exists); 
Ecwl,  the  creeping-wave  field  reradiated  at  point  A;  and 
Ecw2>  the  creeping-wave  field  reradiated  at  point  B. 
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Thus  one  may  write,  the  specular  term  in  the  form 


(10) 


Esp(9)  =  -  1 


COS  0 
j2k0  COS3 


0 

2  - 


JL  e"jk°R 

2R 


and  the  creeping-wave  fields  may  be  written  as 


(ii) 


Ecw  1 


+  Ecw2  =  Ds2  e-j2k°£  [e-(Jk°  +  °s>  0a 

-( jkQ  +  as)  (2tt  -  9)a~l  e_Jk°R 
+  e  -  , 

J  R 


1 1 
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X 

Fig.  4.  Absolute  relative  magnitude  of  backscattered 
fields  due  to  creeping  waves. 


where 


£ 


a  sin 


0 

2 


Using  the  above  equations  the  scattering  cross  section  may  be  written 


(12) 


77  a^ 


1  _ 


:os  0 


Z&S 


■JkO£ 


3  0 

j2k0  cosJ  - 


•(jko  +  as)0a  ,  -(jkQ  +  as)(2Tr- 0)a 
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Fig.  5.  Phase  of  backscattered  fields  due  to  c reeping  waves  . 


A  comparison  of  the  approximate  bistatic  fields  calculated  using 
Eqs  .  (10)  and  (11)  and  the  exact  bistatic  fields  is  presented  in  Fig .  7  . 

It  is  noted  that  these  computations  were  made  by  Moreland,  Peters,  and 
Kilcoyne.  Several  computer  programs  pertaining  to  both  monostatic 
and  bistatic  scattering  by  the  sphere  are  presented  in  Appendix  B.  Results 
obtained  using  these  programs  are  in  agreement  with  the  calculations  of 
Moreland  et  al. 

This  form  of  the  creeping-wave  solution  for  the  sphere  has  the 
advantage  of  a  simple  format.  The  questions  of  the  behavior  of  the 
creeping  wave  fields  in  the  vicinity  of  the  caustic  predicted  by  optics 
at  the  back  of  the  sphere  are  avoided.  Also  it  will  be  seen  that  this 
approach  lends  itself  to  a  straightforward  extension  to  the  analysis  of 
scattering  by  other  bodies  of  revolution.  Furthermore,  the  consistent 
results  obtained  for  monostatic  and  bistatic  radar  cross  section  indicate 
that  there  is  no  need  to  include  an  integration  of  the  creeping-wave  fields 
about  the  shadow  boundary  for  the  back-scatter  case;  i.e.#  to  include  a 
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AT  BlSTATIC  ANGLES,  9,  THERE  ARE  THREE  FIELD 
COMPONENTS  IN  THE  RETURN  SIGNAL.  THEY  ARE 
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Fig.  6.  Bistatic  scattering  of  the  sphere. 


glory  ray  effect.  This  is  consistent  with  the  description  of  the  ray  path 
geometry  given  previously.  The  reader  should  be  cautioned,  however, 
that  this  approach  ensures  success  for  the  back  scatter  case  in  that  the 
effect  of  caustics,  glory  rays,  etc.  can  be  included  in  the  assignment  of 
values  for  the  diffraction  coefficient.  For  this  description  to  be  valid 
requires  that  there  be  no  exceptions  to  it.  A  contradiction  causes  the 
entire  explanation  to  collapse. 


B.  Analysis  of  Backscatter  by  the 
Prolate  Spheroid 


The  prolate  spheroid  is  similar,  with  regard  to  the  form  of  the 
creeping-wave  solution,  to  the  sphere.  As  in  the  case  of  the  sphere,  the 
scattered  fields  are  composed  of  physical  optics  and  creeping-wave  con¬ 
tributions  .  However,  for  the  prolate  spheroid,  the  paths  followed  by  the 
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ig.  7(a).  Real  part  of  the  bistatic  scattered  fields  of  the  sphe 
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Fig.  7(b).  Real  part  of  the  bistatic  scattered  fields  of  the  sphere. 


creeping  waves  on  the  body  are,  in  general,  ellipses.  Thus  the  radius 
of  curvature  in  the  direction  of  propagation  is  continually  changing  and 
it  is  necessary  to  compute  the  integral  in  Eq.  (1).  The  coordinate  system 
used  for  the  prolate  spheroid  is  shown  in  Fig  .  8. 


Fig.  8.  Coordinates  for  the  prolate  spheroid. 


The  prolate  spheroid  has  previously  been  studied  by  Moffatt12  using 
time  domain  analysis.  Moffatt's  approach  is  similar  to  the  creeping-wave 
approach  in  that  he  approximates  the  impulse  response  through  the  use  of 
a  "physical  optics"  contribution  in  the  form  of  an  impulse  at  time  t  =  0, 
which  decays  and  joins  into  a  creeping-wave  return  at  time  t  =  T0  (where 
T0  corresponds  to  the  time  required  for  the  creeping  wave  to  propagate 
around  the  spheroid).  The  form  of  the  impulse  response  postulated  by 
Moffatt  is  shown  in  Fig.  9. 

Following  this  line  of  reasoning  and  referring  to  Fig.  10  we  may 
write  the  back  scattered  creeping-wave  field  for  a  wave  incident  in  the 
x-y  plane,  with  polarization  lying  in  the  x-y  plane,  as 

-\  a(p)ds 

(13)  Ecw  =  2  D2  e"jk°L  e  "A 
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Fig.  9.  Impulse  response  of  the  prolate  spheroid. 


where 


and 

ds  =  differential  arc  length. 

A  similar  expression  may  be  written  for  arbitrary  incidence 
direction  and  polarization,  as  illustrated  in  Fig.  11,  but  we  shall  examine 
the  simplest  case.  For  the  case  illustrated  in  Fig.  10  (parallel  polari¬ 
zation)  we  may  write  the  radius  of  curvature  in  the  incidence  plane  contain¬ 
ing  the  E  vector  as 

(14)  p  =  (a2  sin2  4>  +  b2  cos2  4>)3^2  / ab  , 


which  is  readily  derivable  from  the  parametric  form  of  the  equation  of 
the  ellipse.  Likewise  the  differential  element  of  arc  length  is 
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Using  the  form  of  the  diffraction  and  attenuation  coefficients  given  in 
Eqs  .  (4)  and  (7)  and  leaving  the  magnitude  of  the  coefficients  open  to 
evaluation  using  experimental  data,  one  may  write  the  creeping  wave  as 
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Fig.  11  .  General  incidence  angle  and  polarization 
scattering  by  a  prolate  spheroid. 
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where 


pA  is  the  local  radius  of  curvature  at  point  A  (due  to  symmetry 

Pa=Pb> 

and  the  phase  reference  is  taken  to  be  the  point  (0,  0,  0). 

In  order  to  determine  the  location  of  the  points  of  attachment  and  de¬ 
parture  A  and  B  for  an  arbitrary  incidence  angle  we  must  determine  the 
points  on  the  spheroid  where  the  incident  electric  field  vector  is  parallel 
to  the  normal  to  the  spheroid.  This  is  easily  done  for  the  special  case 
illustrated  in  Fig.  10  through  use  of  the  equation  of  the  surface  and  the 
cross  product  of  the  electric  vector  and  the  normal  vector;  i.e., 


(H) 


and 


1 


a  v 
n  x 


whe  re 


and 


a  .  A  A 

E1  =  -  x  sin  9  +  y  cos  9  . 
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Using  the  above  relations  two  equations  in  two  unknowns  are  obtained, 
the  unknowns  being  the  x  and  y  coordinates  of  the  attachment  points. 

It  is  noted  that  the  attachment  points  are  located  at  (Tx^, +yg)  because 
of  the  symmetry  of  the  body.  Solving,  one  obtains 


(18) 


a2  sin  9 

XA  =  T  Z  Z  I . I 

®  ,Ja2  sin2  0  +  b2  co s 2  0 


and 


(19) 


=  ± 
B 


b2  cos  0 

Ja2  sin2  0  +  b2  cos2  0 


Likewise  it  is  necessary  to  determine  the  point  of  specular  reflection 
so  that  the  proper  phase  between  the  ‘'physical  optics11  contribution  and 
the  creeping  waves  can  be  specified.  This  is  easily  done  using  the 
relations 


(20) 


jsi  +-ii  =  i 


and 

A  A 
v  X  n  =  0, 

where 


A  A 

v  =  x  cos  0 


y  sin  0  ; 


which  results  in 


(21) 


x 


sp 


a2  cos  0 


♦ja2  cos  0  +  b2  sin  0 


and 


(22 )  Ysp  “ 


b2  sin  0 


Ja2  cos  0  +  b2  sin  0 
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Thus,  referring  to  Fig.  10,  the  path  length  difference  between  the 
specular  and  c reeping- wave  components  can  be  written  as 


(23) 


S  =  + 


+  y2 
*  sp 


Lsp 


cos(0-  0Sp)  , 


where 


0Sp  =  arc  tan 


Another  obstacle  to  surmount  in  this  analysis  is  the  fact  that  a 
closed-form  expression  of  the  physical  optics  contribution  for  the 
specular  term  of  the  prolate  spheroid  at  arbitrary  incidence  angles 
does  not  exist.  Moffatt  has  derived  a  time  domain  solution  for  the 
physical  optics  term  by  considering  the  rate  of  change  of  cross  section 
that  a  plane  parallel  to  the  incident  wavefront  encounters  as  it  inter¬ 
cepts  the  spheroid.  This  process  is  illustrated  in  Fig.  12.  Rather 
than  terminating  this  process  at  the  shadow  boundary,  as  is  customary 
for  physical  optics,  Moffatt  terminates  the  physical  optics  integral  when 
the  maximum  cross  section  is  obtained.  Transforming,  Moffatt's  time 
domain  representation  and  discarding  the  contribution  due  to  the  ter¬ 
mination  of  the  integral,  as  is  required  in  creeping-wave  analysis,  one 
may  write  the  "physical  optics"  (i.e.,  the  specular  point)  contribution 


(24) 


Esp 


1  b2  a 

^  b2  sin2  0  +  a2  cos2  0 


j2k>jb2  sin2  0  4*  a2  cos2  0 


•  e 


“  j2k0S 


Computation  of  the  total  field  for  the  prolate  spheroid  is  now  accomplished 
using  Eqs  .  (6)  and  (24).  However,  the  computation  is  not  as  easy  as  for 
that  of  the  sphere  because  of  the  presence  of  the  integrals  in  Eq.  (16). 

In  order  to  perform  the  integration,  numerical  integration  was  employed 
using  a  computer  program  written  to  calculate  the  back  scattered  fields. 
This  computer  program  is  presented  in  Appendix  B.  Also,  it  remains 
to  determine  the  constants  Ul  and  U2  . 
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Fig.  12.  Physical  optics  approximation  for  the  prolate  spheroid. 
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The  computer  program  allows  an  easy  determination  of  the  proper 
values  of  the  parameters  U1  and  U2  through  the  process  of  iteration  and 
comparison. 

It  seems  reasonable  to  expect  that  the  value  Ul  for  the  diffraction 
constant  is  comparable  to  that  of  the  sphere-  As  a  first  approximation 
Ul  is  taken  equal  to  the  value  for  the  sphere,  i-e.,  (0.27).  Also  it 
seems  reasonable  that  the  value  for  the  attenuation  coefficient  should 
be  reasonably  close  to  that  of  the  sphere  for  spheroids  of  moderate 
eccentricities.  The  case  of  the  2:1  spheroid  is  examined  here  since 
a  large  amount  of  data  are  available  for  this  target.  It  is  of  interest 
to  compare  the  maximum  and  minimum  radii  of  curvature  for  a  2:1 
prolate  spheroid.  It  is  apparent  that  the  maximum  and  minimum  radii 
of  curvature  in  the  x-y  plane  occur  at  q>  =  0  and  tt/2,  respectively. 

Using  Eq .  (14)  we  find  that 


(25)  pmin  ~ 

and 


8 

If  we  define  the  ratio  of  the  radius  of  curvature  in  the  direction  of 
propagation  of  the  creeping  wave  to  the  orthogonal  radius  of  curvature 
to  be  Rq,  we  have  the  following  correspondence: 

Body  Ra  U2 

Cylinder  0  1.68 

Sphere  1  0.84 


p  —  4b  , 

Kmax 


thus  the  ratio  pmin/  pmax  = 


2:1  Prolate  Spheroid 
(Parallel  Polarization) 


^■Ohnin  ~ 


Not  yet 
determined 


Thus  one  expects  that  the  value  of  U2  for  the  prolate  spheroid  will  be 
equal  to  or  less  than  that  of  the  sphere,  and  may  indeed  vary  along  the 
surface  of  the  spheroid.  However,  to  a  first  approximation  a  constant 
value  of  U2  will  be  assumed,  since  the  form  of  a  possible  correction 
term  to  the  attenuation  coefficient  is  not  known. 
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Using  the  computer  program,  the  back  scatter  as  a  function  of 
incidence  angle  was  computed  for  a  range  of  U2  values  for  0.2  to  0.8, 
and  with  Ul  equal  to  0.27.  The  results  of  these  computations  are  shown 
in  Fig.  13.  It  is  seen  that  for  kQa  values  greater  than  approximately 
3.0,  satisfactory  results  are  obtained.  However,  the  present  form  of 
the  analysis  results  in  a  creeping  wave  contribution  which  is  too  large 
for  values  of  kQa  less  than  3.0.  As  a  further  check  computations  over 
a  range  of  koa  values  were  made  for  the  special  case  (perpendicular 
polarization)  shown  in  Fig.  14.  Here  we  see  that  the  path  traversed 
by  the  creeping  wave  is  circular  with  a  value  Ra  =  1/2,  and  U2  was 
chosen  to  be  0.84  as  a  first  approximation.  An  adaptation  of  the  previous 
computer  program  for  the  sphere  was  made  by  using  the  ’’physical  optics" 
form  for  the  spheroid.  The  resultant  program  is  presented  in  Appendix  B. 
Figure  15  presents  the  results  of  this  computation,  and  it  is  seen  that 
excellent  agreement  between  the  computed  and  measured  values  is  obtained. 
It  is  also  noted  that  the  curve  computed  using  creeping-wave  theory  is  the 
same  as  (within  computational  accuracy)  the  curve  predicted  by  Moffatt 
using  time  domain  analysis. 
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Fig.  13.  Backscatter  patterns  of  the  prolate  spheroid  for  a 
range  of  U2  values  (parallel  polarization). 
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Fig.  13.  Backscatter  patterns  of  the  prolate  spheroid  for  a 
range  of  U2  values  (parallel  polarization). 
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Fig.  14.  Scattering  by  the  prolate  spheroid  for  perpendicular 
polarization  and  0  =  90°  . 

From  the  excellent  results  of  Fig.  15  it  would  seem  that  a  more 
accurate  calculation  of  the  patterns  of  Fig.  13  should  be  possible  for 
small  values  of  kQa .  The  inaccuracies  in  Fig.  13  for  small  values  of 
kQa  are  seen  to  occur  in  the  region  about  0  =0  (i.e.,  nose -on).  A 
possible  explanation  of  these  inaccuracies  is  that  the  effects  of  the  small 
radius  of  curvature  at  the  "nose"  of  the  spheroid  are  not  accurately  pre¬ 
dicted  by  a  single  term  in  the  attenuation  coefficient.  Another  possibility 
is  that  a  reflection  occurs  at  or  near  the  nose  thus  decreasing  the  effect 
of  the  creeping-wave  fields.  Future  work  should  be  devoted  to  extension 
of  the  analysis  to  smaller  spheroids  and  the  formulation  of  a  computer 
program  for  a  general  angle  of  incidence.  A  goal  of  this  research  would 
be  to  evaluate  the  constant  associated  with  the  attenuation  factor  as  a 
function  of  the  ratio  of  the  principal  radii  of  curvature;  i.e.,  a  more 
accurate  form  for  U2  vs.  Ret  should  be  found.  The  values  for  Ra  =  0, 

1,  and  oo  have  already  been  established  within  reasonable  limits. 


28 


2: 1  PROLATE  SPHEROID 


r 


T  lyr'^TW”*'  rr][rr 


(2\  < 9P)  -o 


29 


Fig.  15a.  Echo  area  of  the  prolate  spheroid  as  a  function  of  k, 
for  perpendicular  polarization  and  9  =  90°  . 
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for  perpendicular  polarization  and  0  =  90 


C.  Analysis  of  the  Ogive  and  Spherically 
Capped  Ogive  in  the  Near  Nose-On 
Region 


The  backscatter  of  the  ogive  has  been  treated  by  Peters1  using  a 
traveling  wave  antenna  approach.  A  modification  of  this  traveling  wave 
analysis  has  been  developed  by  Peters11  which  accounts  for  radiation 
losses  through  use  of  a  creeping-wave  attenuation.  Cohen7  experiment¬ 
ally  determined  values  for  the  reflection  at  a  tip  of  the  ogive.  Here  we 
shall  show  that  a  consistent  picture  of  scattering  by  the  ogive  can  be 
formulated  using  creeping-wave  analysis  and  application  of  geometrical 
theory  of  diffraction  techniques.  In  particular  it  is  shown  that  the  wedge 
diffraction  coefficient  derived  by  Pauli18  is  applicable  to  the  determina¬ 
tion  of  the  tip  reflection  and  transmission  coefficients. 

Figure  16  shows  the  ogive  and  the  coordinate  system.  As  in  the 
case  of  the  prolate  spheroid,  incident  energy  is  trapped  at  points  A  and 
B  where  the  incident  electric  vector  is  parallel  to  the  normal  to  the 
surface.  These  creeping  waves  travel  along  the  ogive  until  they  en¬ 
counter  the  tip. 

At  the  tip,  part  of  the  creeping- wave  energy  is  diffracted  as 
creeping  waves  on  all  geodesics  containing  the  tip  and  the  remainder 
is  radiated  from  the  tip#  The  previous  treatment11  assumed  that  the 
creeping-wave  fields  on  all  of  the  geodesics  have  equal  magnitudes. 

This  is  approximately  correct  but  the  small  differences  in  the  magni¬ 
tudes  will  yield  the  on-axis  c reeping-wave  field  which  can  be  observed 
in  the  measurements  reported  by  Blore.19  Off-axis  the  creeping  wave 
return  is  a  result  of  the  creeping  waves  which  propagate  on  the  ogive 
in  the  plane  defined  by  the  incidence  direction  and  the  incident  E-field. 

In  this  plane  the  effect  of  the  tip  can  be  characterized  by  reflection  and 
transmission  coefficients  at  the  tip.  Siegel  et  al  have  derived  an  ex¬ 
pression  for  the  physical  optics  scattering  by  the  ogive  by  assuming  that 
for  both  monostatic  and  bistatic  returns  in  the  near  nose-on  region  the 
illuminated  tip  is  the  dominant  scatterer.  Their  expression  for  on-axis 
backscatter  is 


(26) 
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Referring  to  Fig.  16,  the  backscatte red  fields  attributed  to  the 
creeping  waves  may  be  written  as 
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where 

D  is  the  ogive  diffraction  coefficient, 

Rt  is  the  tip  reflection  coefficient, 

Tt  is  the  tip  transmission  coeffection, 
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aQg  is  the  ogive  attenuation  coefficient,  and 
S,  L1}  L2  are  as  shown  in  Fig.  16. 

Bistatic  radar  cross  section  measurements  of  the  ogive  have  shown 
that  the  energy  incident  on  the  tip  is  scattered  by  the  tip  along  each 
geodesic  on  the  ogive;  i.e.,  these  tip- scatte red  fields  show  no  4> 
dependence.  This  agrees  with  the  results  of  the  equivalent  antenna 
approach,  for  which  no  4>  variation  of  the  field  was  assumed.  Thus 
a  factor  of  1/n/2tt  is  introduced  to  account  for  the  lack  of  j>  dependence 
of  the  tip- scatte red  fields.  Again  attachment  points  A  and  B  are 
determined  by  finding  the  points  on  the  surface  where  the  normal  to 
the  surface  is  parallel  to  the  incident  electric  vector.  Straight¬ 
forward  analysis  results  in  the  relations 


(28)  S  =  h  sin  0, 

Ll  =  r  [  _£  .  0  ]  (  and 


L2  =  r[_?E.  +  0 


In  the  case  of  the  ogive  Peters11  has  postulated  that  the  diffraction 
coefficient  remains  unchanged  from  that  of  the  sphere  and  that  the 
attenuation  coefficient  is 


TT 

(29)  aog  =  i  as  =  0.21  \S  e  6 

It  is  proposed  here  that  the  reflection  and  transmission  coefficients 
be  determined  using  wedge  diffraction  theory  and  geometric  diffraction 
theory . 

It  is  seen  from  Fig.  17a  that  the  tangent  planes  to  the  surface  at  the 
tip  of  the  ogive,  which  are  normal  to  the  plane  defined  by  the  incidence 
direction  and  the  electric  field  vector,  form  a  wedge.  Using  the  wedge 
diffraction  formulation  of  Pauli18  we  write 


(30) 


Rt  =  n/2tt  kr  e^r  vg(r,  n,  rr) 
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Fig.  17  .  Tangent  planes  at  the  tip  of  the  ogive  . 


and 

Tt  =  \l2w  kr  e^r  vg(r,  n,  mr)  , 

where 


n  =  —  (2tt  -  WA) 

TT 

and  where  the  form  of  vg  given  in  Eq.  (48)  for  large  values  of  kQr  is 
employed.  Only  the  relative  magnitude  of  the  diffracted  fields  is  utilized, 
resulting  in  the  multiplication  of  vg  by  the  factor  n/2tt kr  e .  Previously, 
Cohen7  had  used  an  experimental  approach  to  determine  values  of  the 
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reflection  coefficient.  His  method  actually  determined  the  sum  of  the 
reflected  and  transmitted  components. 

The  radar  cross  section  for  long  thin  bodies  has  been  derived  by 
considering  the  body  to  be  an  equivalent  lossless  antenna.1  The  radar 
cross  section  of  an  antenna  is 


(31) 
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where 


G(9,0>)  is  the  gain  of  the  antenna, 

0, 9  are  .he  far  field  angles,  and 

y  is  the  voltage  reflection  coefficient  at  the  distant  end. 

The  approximate  echo  area  obtained  by  this  equivalent  antenna  approach 
is  compared  in  Fig.  18  to  the  experimental  echo  area.  The  difference 
between  this  lossless  theoretical  picture  and  experimental  values  is 
attributed  to  a  loss  mechanism.  This  loss  mechanism  could  not  result 
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from  energy  absorption  since  the  target  is  nearly  perfectly  conducting. 
It  has  been  shown11  that  this  loss  mechanism  is  due  to  radiation,  and 
that  the  attenuation  coefficient  appropriate  to  determining  the  loss  is 
that  of  Eq.  (29).  Thus  the  geometrical  theory  of  diffraction  solution 
can  be  used  in  the  previous  equivalent  antenna  approach  to  correct  the 
lossless  radar  cross  section  computation.  Computing  the  gain  pattern, 
G(Q,4>)>  as  in  the  previous  analysis11*7  and  using  the  ogive  attenuation 
coefficient  to  correct  the  gain  for  radiation  loss,  one  obtains 


(32) 


°los  s 
°"lossless 


The  resultant  cross  section  computed  using  this  correction  for  radiation 
loss  is  shown  in  Fig.  18  and  is  seen  to  be  in  excellent  agreement  with 
the  experimental  results  . 

Equation  (27)  has  been  used  to  calculate  the  maximum  scattered 
fields  for  ogives  having  tip  angles  of  40°  and  30°  with  radii  of  1  1 . 3X 
and  77\,  respectively.  A  result  is  shown  in  Fig.  18.  It  is  seen  that 
satisfactory  agreement  was  obtained  in  the  case  of  the  40°  ogive  but 
the  agreement  in  the  case  of  the  30°  angle  was  not  as  good.  In  addition, 
the  scattered  fields  on-axis  were  computed  using  Eq.  (33); 


(33) 
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Equation  (33)  for  the  on-axis  field  was  obtained  by  assuming  that  the 
form  of  Eq.  (27)  holds  and  that  the  fields  diffracted  by  the  tip  in  the 
shadow  follow  geodesics  on  the  ogive.  Also  it  is  assumed  that  the  tip- 
scattered  fields  are  uniform  thus  producing  a  radial  field  distribution 
at  the  shadow  boundary.  Integration  of  the  contributions  around  the 
shadow  boundary  produces  the  factor  of  two  difference  in  Eq.  (27)  and 
(33).  Results  obtained  using  Eq.  (33)  were  compared  with  the  results 
obtained  by  Blore19  for  an  ogive  with  a  tip  angle  of  40°  and  a  radius  of 
5.65X.  It  was  found  that  Eq.  (33)  did  not  accurately  predict  the  axial 
creeping- wave  contribution  because  the  magnitude  of  |t|  -  R  was 
too  large.  The  results  of  Blore  indicate  that  the  on-axis  creeping- wave 
contribution  must  be  of  the  same  order  of  magnitude  as  the  tip- scattered 
field,  which  may  be  expressed  as 
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(34) 


The  above  results  indicate  that  Rt  and  Tt  are  nearly  equal,  a  fact 
not  predicted  by  wedge  diffraction. 

Unfortunately  Rt  and  Tt  are  evaluated  for  two-dimensional  geometry. 
There  is  no  canonical  solution  for  these  diffraction  parameters  when 
the  wedge  is  replaced  by  a  cone  tip  canonical  surface.  However  this 
approach  does  explain  the  presence  of  the  creeping-wave  term  for  the 
ogive  on-axis  in  that  Rt  and  Tt  are  different  in  magnitude,  thus  result¬ 
ing  in  an  on-axis  creeping-wave  contribution,  although  the  resulting 
contribution  is  in  error.  Further  theoretical  and  experimental  exami¬ 
nation  of  this  problem  is  required  in  order  to  obtain  a  more  accurate 
solution  for  the  ogive. 

Moreland,  Peters,  and  Kilcoyne  have  applied  c reeping- wave  theory 
to  analysis  of  a  spherically  capped  ogive  as  shown  in  Fig.  19-  They  have 
expressed  the  scattered  field  as 


(35) 
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=  radius  of  curvature  of  the  ogive, 

=  radius  of  the  spherical  cap, 

=  angles  shown  in  Fig.  19, 

=  attenuation  coefficient  of  ogive, 

=  diffraction  coefficient, 

=  Jko  +  as  ’  and 

=  attenuation  coefficient  of  the  sphere. 
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The  reflection  coefficient,  T,  and  the  transmission  coefficients,  T0s 
and  Tso>  arise  because  of  the  change  of  surface  impedance  at  the  junction 
of  the  ogive  and  sphere.  The  surface  impedance  of  the  ogive  is  taken  to 
be  ZQ,  the  impedance  of  free  space,  and  the  surface  impedance  of  the 
sphere  is  taken  to  be  the  ratio  of  E0  to  HQ  for  the  fields  of  the  major 
creeping  waves  postulated  in  Eqs  .  (5)  and  (6)  for  the  sphere.  Referring 
to  Fig.  3  we  see  that  this  ratio  is  dependent  upon  the  sphere  radius  and 
can  be  determined  from  the  exact  fields  of  the  sphere  for  a  given  radius. 
Thus  T  and  Tos  can  be  written  as 


(36)  r  =  Zg  ~  Z° 

Z s  ^  Zo 


and 
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Moreland  et  al8  found  that  it  was  necessary  to  assign  a  negative  value  to 
T  in  order  to  obtain  meaningful  results,  but  offered  no  explanation  for 
this  phenomenon.  The  explanation  for  the  difference  in  sign  between  the 
terms  in  Eq  .  (35)  is  illustrated  in  Fig.  20  and  is  explained  using  geo¬ 
metrical  theory  of  diffraction  analysis.  It  is  known  that  if  a  ray  path 
encounters  a  discontinuity  in  a  surface  (as  shown  in  Fig.  20  )  the 
incident  ray  is  diffracted  at  the  discontinuity  giving  rise  to  the  exist¬ 
ence  of  a  caustic  at  that  point.  The  case  of  a  "ring"  discontinuity,  as 
exists  here,  introduces  a  phase  shift  of  tt/  2  radians  in  the  ray  passing 
through  the  discontinuity.  Figure  20  ,  due  to  Kouyoumjian,  15  illustrates 
this  process.  If  the  process  of  reflection  and  transmission  of  the  incident 
ray  is  described  as  follows,  it  is  seen  that  this  phase  shift  must  be  taken 
into  account.  The  process  may  be  described  in  the  following  steps-  The 
incident  ray,  which  is  bound  to  the  surface,  encounters  the  discontinuity 
and  gives  rise  to  reflected  and  radiated  fields.  The  reflected  field  is 
obtained  using  Eq .  (36).  Next  the  radiated  field  in  the  direction  of  the 
surface  is  obtained  from  the  transmission  coefficient.  The  transmitted 


Fig.  20.  Diffraction  at  a  surface  discontinuity 
which  gives  rise  to  a  caustic  . 
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field  proceeds  through  the  discontinuity  undergoing  a  phase  shift  due  to 
the  caustic  of  the  radiated  rays  and  re-attaches  itself  to  the  surface. 
Since  the  transmitted  field  undergoes  this  process  twice  in  its  travel 
around  the  spherically  capped  ogive,  the  total  phase  shift  is  tt  radians. 

In  addition,  the  factor  1/n/Ztt  proposed  by  Moreland  et  al  is  retained  to 
account  for  the  rotational  symmetry  of  the  caustic  which  causes  dif¬ 
fracted  rays  to  propagate  along  the  spherical  cap  in  directions  other 
than  that  of  the  plane  of  incidence,  thus  decreasing  the  backscattered 
field.  Combining  all  the  terms  the  total  expression  for  the  backscat¬ 
tered  field  is 

(38)  Es  =  Ei  Dog  e"2Yogr0o  jr 

l 

-j(k0«  +  Ysrs  2  0i+  j2tt) 

The  results  of  the  computation  for  the  spherically  capped  ogive  are 
presented  in  Fig.  21  and  compared  to  experimental  measurement. 

It  is  seen  that  in  the  near  nose-on  region  the  results  are  in  good  agree¬ 
ment  with  measurement.  This  form  of  creeping-wave  analysis,  coupled 
with  the  traveling  wave  picture  and  physical  optics,  allows  prediction  of 
the  backscattered  field  of  the  ogive  over  a  wide  range  of  incidence 
directions  and  k0a  values  .  It  should  be  noted  that  the  impedance  of  the 
creeping  wave  is  obtained  by  a  best  fit  of  the  fields  at  the  surface  of  the 
sphere  in  the  entire  shadow  region.  This  would  include  the  fields  close 
to  the  point  at  the  rear.  It  has  been  noted  that  this  is  a  region  of  poor 
fit,  regardless  of  the  approximation  used,  because  there  are  fields  pre¬ 
sent  here  which  are  not  included  in  the  simple  creeping-wave  picture. 
This  creeping  wave  picture,  as  noted  previously,  does  not  include  a 
caustic  at  the  point  at  the  rear.  Further  improvement  in  these  results 
would  be  anticipated  if  the  fields  were  more  closely  matched  by  two 
creeping  waves  at  a  distance  removed  from  the  rear  point  and  if  multiple 
interactions  are  included. 
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II.  CREEPING  WAVE  ANALYSIS  OF  PLANAR 
TARGETS  FOR  EDGE- ON  INCIDENCE 

A.  The  Circular  Disk 

Ufimtsev^0  and  DeVore  and  Kouyoumjian21  have  treated  the  circular 
disk  for  arbitrary  incidence  angle  and  polarization  using  diffraction  tech¬ 
niques.  For  the  case  of  edge-on  incidence  with  polarization  parallel  to  the 
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A  -  CAP  RADIUS  ■  0.007“,  90  •  4*{  A  ^ppRQX. 
0-  CAP  RADIUS  •  0.695",  90  »  8°  {” 0"”  APPROX 


ANGLE  OF  INCIDENCE 
(a) 

C-  CAP  RADIUS-  0.509“  90  •  '2*  ^ ” APPROX. 
D  -  CAP  RADIUS  •  0.254"  90  •  16°  approx 


AN6LE  OF  INCIDENCE 

(b) 


Fig.  21  .  Backscattered  field  pattern  of  the 
spherically  capped  ogive. 
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edge  (as  shown  in  Fig.  22)  it  is  possible  to  construct  a  solution  using 
creeping-wave  concepts.  First  the  specular  return  is  computed  by 
applying  geometrical  optics  in  the  plane  of  the  disk  and  edge  diffraction 
in  the  plane  normal  to  the  edge.  As  seen  from  Fig.  22,  the  geometrical 
optics  spreading  factor  in  the  plane  of  the  disk  can  be  found  by  considera¬ 
tion  of  the  reflected  rays.  It  is  seen  that  these  reflected  rays  diverge 
from  the  virtual  focus  located  a  distance  a/  2  from  the  center  of  the  disk. 
The  geometrical  optics  spreading  factor  is  written  as 

[( P 1  P2/(Pl+  *)(P2  +  0]  2.  15 


where  for  plane  wave  incidence  •  In  this  case  p2  is  the  focal  i 

distance  a/ 2  and  thus  the  space  attenuation  factor  becomes  (a/(a+2£))2. 
Combining  the  space  attenuation  factor  and  the  result  using  the  edge 
diffraction  coefficient  of  Appendix  A  we  have,  for  the  backscattered 
field  from  the  specular  point, 


Fig.  22.  Scattering  by  a  circular  disk  edge  on. 
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Next  it  is  postulated  that  two  creeping  waves  travel  around  the  disk,  as 
in  the  case  of  the  cylinder  and  sphere.  Thus  the  creeping-wave  back- 
scattered  field  can  be  written  in  the  form 


(40)  Ecw  =  -  213^  e'jZk°a  e'(jk°  +  Qd)lTa 


It  remains  to  determine  the  attenuation  and  diffraction  coefficients  which 
apply  to  such  a  target.  If  the  disk  is  considered  to  be  the  limiting  case  of 
an  oblate  spheroid,  it  is  seen  that  the  ratio  of  the  orthogonal  radii  of 
curvature  at  the  edge  of  the  disk  can  be  expressed  as 


(41)  Ra  =  ^ 

b 


where  a  and  b  are  shown  in  Fig.  23. 

Now,  as  the  oblate  spheroid  is  "squashed  down"  into  a  disk  and 
b  —  0  it  is  seen  that  the  ratio  R  —+oo .  Relating  this  to  the  previous  ratios 
of  the  cylinder,  sphere,  and  ogive  it  is  reasonable  to  expect  that  the 
attenuation  coefficient  will  possess  a  magnitude  U2  on  the  order  of  or 
smaller  than  that  of  the  ogive.  As  a  first  trial  a  value  U2  =  0.21  (i.e., 
that  of  the  ogive)  was  chosen.  It  is  more  difficult  to  speculate  about  the 
form  and  magnitude  of  so  the  experimental  data  of  DeVore  and 
Kouyoumjiam  have  been  used  to  determine  D^.  The  experimental  data 
indicate  a  deep  null  in  the  cross  section  for  a  value  of  kQa  =  6.  Using 
Eq.  (4u)  and  taking  =  aQg  the  total  field  was  set  equal  to  zero  for 
kQa  =  6.  Next  it  is  postulated  that  D^2  is  of  the  same  functional  form 
as  previously  stated.  This  calculation  results  in 

L  2  -j  _JL  +  •  IjL  -i-  •  — 

(42)  Dd2  =  0.0974  a-K-S  e  12  e  J  4  ^  j  D2phere  e  J  4  . 


This  value  of  D2  was  then  used  in  connection  with  the  computer  program 
in  Appendix  B  to  calculate  the  curve  shown  in  Fig.  24.  In  order  to  check 
the  results  for  kQa  values  larger  than  those  reported  by  DeVore  and 
Kouyoumjiam,  additional  measurements  were  made.  These  measurements 
are  shown  as  a  range  of  values  because  the  difficulty  in  obtaining  great 
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Fig,  23.  The  oblate  spheroid  as  an  approximation  to  the  disk. 


precision  in  the  measurements  of  such  low  cross  section  targets.  It  is 
seen  that  this  analysis  yields  results  which  are  in  surprising  agreement 
with  measurement.  It  seems  reasonable  to  examine  other  edge-on  con¬ 
figurations  to  determine  if  significent  creeping-wave  contributions  are 
caused  by  edges. 


B.  Analysis  of  the  Flat  MOgivel! 


The  analysis  of  the  flat  ogive  parallels  that  of  the  ogive  of  revolution 
with  the  attenuation  and  diffraction  coefficients  developed  for  the  disk  sub¬ 
stituted  for  the  previous  ogive  attenuation  and  diffraction  coefficients. 

Again,  Pauli's  form  for  wedge  diffraction  is  employed  to  specify  reflection 
and  transmission  coefficients  at  the  tip.  However,  for  the  case  of  the  flat 
ogive  the  tip  is  treated  as  a  wedge  rather  than  as  a  point  in  the  determination 
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Fig-  24.  Backscattering  by  the  circular  disk  edge  one 
for  a  range  of  kQa  values  ♦ 
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of  the  geometrical  optics  spreading  factor.  The  wedge  diffraction  coef¬ 
ficient  given  in  Appendix  A  accounts  for  this  spreading  of  rays.  The 
resulting  expression  for  the  scattered  field  is 


(43)  Es  =  El  Dd 


o  +  aog)  e-j2k0s 


d-  e-(jko  4  ^og)^^-12 


[' 


e-(jkG  +  aog)(L!  + 


e 


-jkR 


where 


Dd  is  given  by  Eq.  (42)  , 


R\y  and  T\y  are  given  by  Eq.  (30)  , 


aQg  is  the  ogive  attenuation  coefficient,  and 
L»l,  Lz,  and  s  are  shown  in  Fig.  16  . 


Equation  (43)  has  been  used  to  calculate  the  echo  area  of  the  flat  ogive 
and  the  results  are  presented  in  Fig.  25.  It  is  noted  that  the  maximum 
return  of  the  flat  ogive,  as  measured  and  calculated,  is  greater  than  the 
return  of  the  ogive  of  revolution.  This  is  a  result  of  the  difference  in 
the  geometrical  optics  spreading  at  the  tip.  Whereas  the  ogive  tip  dif¬ 
fracts  energy  into  all  sectors  of  three-dimensional  space,  the  flat  ogive 
wedge  concentrates  the  diffracted  energy  in  the  plane  of  the  flat  ogive. 
Thus  the  geometrical  optics  spreading  factors  at  the  tip  differ  by  Ztt 
which  compensates  for  the  lower  diffraction  coefficient  of  the  flat  ogive 
and  results  in  a  higher  return  for  the  flat  ogive.  This  rather  surprising 
fact  is  substantiated  by  the  measurements  shown  in  Fig.  25  and  it  is  seen 
that  good  agreement  between  theory  and  experiment  is  obtained. 

The  results  of  the  analysis  of  the  flat  ogive  indicate  that  large 
creeping-wave  contributions  to  the  scattered  field  of  a  target  can  result 
from  the  presence  of  edges.  This  type  of  creeping- wave  contribution  is 
limited  to  the  plane  containing  the  "fin11.  In  scattering  directions  where 
the  creeping  wave  on  the  edge  does  not  contribute,  it  is  possible  to  obtain 
an  estimate  of  the  scattered  field  using  wedge  diffraction  techniques.22*23 
At  the  present  time  only  an  estimate  is  possible  because  the  diffraction 
coefficient  for  a  finite  edge  is  not  known.  This  estimate  is  obtained  by 
applying  the  diffraction  coefficient  in  the  plane  normal  to  the  edge  and 
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Fig.  25.  Backscatter  of  the  ,!flatM  ogive. 
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assuming  that  the  edge  acts  as  a  traveling  wave  antenna  with  constant 
current  excitation  in  the  plane(s)  containing  the  edge.23  Such  an 
assumption  does  not  meet  the  boundary  condition  on  a  finite  edge  but 
does  provide  an  estimate  of  the  maximum  scattered  field. 

Further  study  of  the  effects  of  creeping  waves  on  edges  in  com¬ 
bination  with  boundary  value  techniques  may  lead  to  a  diffraction  type 
solution  for  scattering  by  a  finite  edge.  Such  a  solution  will  allow 
rapid  calculation  of  the  effects  of  the  presence  of  fins  and  edges  on  a 
target . 


III.  CONCLUSIONS 

It  has  been  demonstrated  that  the  geometrical  theory  of  diffraction 
which  incorporates  edge  diffraction,  creeping  wave,  and  geometrical 
optics  techniques  can  be  applied  to  the  problem  of  scattering  by  volu¬ 
metric  bodies  and  edges.  In  the  sample  cases  considered  (spheres, 
prolate  spheroids,  ogives,  capped  ogives,  disks,  and  flat  ogives) 
solutions  for  the  backscatte red  fields  were  obtained  in  simple  analytical 
form.  Computation  of  the  backscatte  red  fields  using  these  simple  forms 
is  straightforward  and  takes  but  a  small  amount  of  time  on  a  computer. 

The  nature  of  the  geometrical  theory  of  diffraction  solution  makes  it 
easy  to  determine  the  major  scattering  centers  on  a  target  through  the 
use  of  ray-tracing  techniques  and  simple  computations.  Thus  a  rapid 
estimate  of  radar  cross  section  can  be  obtained  for  a  special  case. 

Current  limitations  on  this  method  include  the  lack  of  a  diffraction 
solution  for  scattering  by  a  finite  edge,  and  the  unknown  dependence  of 
the  attenuation  and  diffraction  coefficients  on  the  orthogonal  radius  of 
curvature  of  a  solid  body.  However,  results  obtained  to  date  indicate 
that  the  effects  of  the  orthogonal  radius  of  curvature  may  be  approximated 
with  good  accuracy.  It  also  appears  that  no  problem  exists  for  which 
continued  research  would  not  yield  a  solution.  A  satisfactory  approxi¬ 
mation  for  the  scattering  by  a  finite  edge  seems  possible  using  a  combina¬ 
tion  of  diffraction  and  boundary  value  techniques.  Any  problem  that  cannot 
be  solved  analytically  can  be  solved  by  experimental  techniques;  e.g.,  the 
reflection  and  transmission  coefficients  for  the  ogive.  The  unifying  char¬ 
acteristics  of  the  method  offer  the  expectation  that  the  shadow  region  geo¬ 
metry  of  more  complex  targets  could  be  included  in  radar  cross  section 
solutions,  or  even  that  the  characteristics  of  the  shadow  region  could  be 
determined  from  the  measured  scattered  fields  caused  by  the  electro¬ 
magnetic  energy  that  has  propagated  through  this  region. 
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IV.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


The  concept  of  creeping  wave  diffraction  has  been  established  in 
this  report  for  geometries  that  have  not  been  previously  considered 
from  this  point  of  view.  It  remains  to  obtain  more  precise  values  of 
the  diffraction  coefficients  and  attenuation  factors  for  the  shapes  dis¬ 
cussed  herein.  These  efforts  should  not  be  restricted  to  any  one  ap¬ 
proach.  Experimental  data  and  modern  computational  methods  such 
as  the  point  matching  boundary  value  techniques  as  well  as  the  concepts 
of  wedge  diffraction  should  be  used  to  obtain  improved  values  for  these 
quantities  . 

Applications  of  wedge  diffraction  and  the  impedance  of  the  creep¬ 
ing  wave  should  be  made  to  evaluate  the  effects  of  a  discontinuity  in  the 
shadow  region.  Attention  should  be  focused  on  the  exact  fields  at  the 
surface  of  the  sphere.  A  serious  attempt  should  be  made  to  reconstruct 
these  fields  from  the  creeping  wave  point  of  view.  Such  studies  should 
yield  a  more  profound  understanding  of  the  creeping  wave  mechanisms 
and  the  limitations  of  this  approach.  The  study  should  be  extended  to 
geometries  other  than  those  discussed  in  this  report.  An  ultimate  goal 
of  this  approach  would  be  to  reduce  the  computation  of  the  scattered 
fields  of  this  type  to  the  level  that  practicing  engineers  could  readily 
evaluate  them  or  conversely  that  a  generalized  computer  program  could 
perform  this  task  for  a  body  of  general  shape. 

The  applications  to  antenna  theory  should  not  be  overlooked.  An 
aircraft  is  a  curved  surface.  An  antenna  on  such  a  surface  could  be 
effected  by  such  a  creeping  wave.  Ground  planes  could  be  terminated  in 
a  curved  edge.  This  would  reduce  irregularities  that  appear  in  antenna 
patterns  mounted  on  such  a  ground  plane.  In  this  case  it  becomes 
necessary  to  separate  the  diffraction  coefficient  into  its  two  components 
as  trapping  and  launching  mechanisms  do  not  necessarily  occur  together, 
i.e.,  the  antenna  parameters  would  replace  one  or  the  other  of  these  on 
transmission  or  reception  respectively. 
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APPENDIX  A 


In  order  to  satisfactorily  calculate  diffracted  rays  it  is  necessary 
to  determine  the  diffraction  of  plane  and  cylindrical  waves  by  a  wedge* 

This  has  previously  been  done  (see  Ref.  24)  and  an  expanded  treatment 
is  presented  here.  Sommerfeld  obtained  the  solution  for  a  perfectly  con¬ 
ducting  wedge  composed  of  two  half-planes  with  a  plane  wave  incident  on 
the  wedge.  Sommerfeld  also  obtained  an  explicit  form  of  the  solution  for 
zero  wedge  angle;  i.e.,  a  half-plane.  Pauli  determined  an  explicit 
form  for  the  general  wedge  .  Oberhettinger25  presents  a  different  form 
for  the  general  wedge  which  he  obtained  using  Green's  function  techniques. 

The  diffraction  of  a  cylindrical  wave  by  a  wedge  is  obtained  by  the 
use  of  reciprocity,  together  with  Pauli's  expressions  for  wedge  diffraction. 
Pauli’s  expressions  give  the  diffraction  of  a  plane  wave  by  a  wedge  for  a 
general  angle  of  incidence  and  for  polarization  either  perpendicular  or 
parallel  to  the  edge  of  the  wedge.  The  total  field  at  observation  point  P 
of  cylindrical  co-ordinates  (r,i)j),  as  shown  in  Fig.  26,  is  given  by 


Fig.  26.  Diffraction  by  a  wedge  of  included  angle  (2-n)iT. 
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(44) 


u(r,  n,  i|i)  =  v(r,  n,  ili-vp0)  ±  v(r,  n>  4*  +4'0) 


for  a  plane  wave  incident  from  direction  i|/Q.  The  plus  sign  applies  for 
the  polarization  of  the  electric  field  perpendicular  to  the  edge 


and  the  minus  sign  applies  for  polarization  parallel  to  the  edge 
The  quantities  v(r,4>)  are  given  by 


=  0 


)• 


(45)  v(r,  n,  9)  =  v*(r,  n,  4>)  +  vs(r,  n,  9)  , 


where  vv ( r ,  n,  $)  is  the  geometrical  optics  field;  and  is  given  by 


(46) 


where 


vv(r,  n,  cj>)  = 


j  exp[  jkr  cos  (cj)  +  2tt  nN)]  ,  -tt  +  2tt  nN  <  (j>  <  tt  +  2tt  nN 
0  otherwise, 


i 


N  =  0,  +1,  ±2,  ... 


vjj(r,  n,  4>)  is  the  diffracted  field  for  a  wedge  of  angle  WA 
given  by 


(47) 


vB(r,  n,  9) 


2e 


TT 

J4 


'Vrr 


4 

cos  ~ 

1 

2 

TT  q) 

COS  —  COS 


(2-ii)tt;  and  is 


^jkr  cos  4> 


dT 


where  a  =  1  +  cos  $  . 
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Equation  (47)  is  composed  of  a  leading  term  plus  higher-order 
terms  which  are  negligible  for  large  values  of  kr.  For  large  values  of 
(akr)  Eq.  (47)  becomes 


1  t r 

—  s  m  — 

n  n  . 

■  —  - _ — _ _  ,  T  «  •  «  • 

T T  <j> 

COS  —  -  COS  — 

n  n 

The  diffracted  field,  as  expressed  by  Eq.  (48),  is  that  from  which  the 
asymptotic  diffraction  coefficients  of  the  geometrical  theory  of  dif¬ 
fraction  are  obtained.  Since  this  expression  is  valid  only  for  large 
values  of  (akr),  it  is  not  valid  on  the  shadow  boundary  because  a  =  0 
there.  Then  Eq .  (47)  must  be  used,  which  gives  the  value  of  the  dif¬ 
fracted  field  on  the  shadow  boundary  as 

(49)  Vg(r,  n>  ir)  ~  ^  J  e  ^ r  +  .  .  . 

upper  sign  for  <j>  =  tt  “ 
lower  sign  for  <j>  =  t r+ 

The  value  of  v(r,c)))  on  the  shadow  boundary  can  then  be  obtained 
from  Eqs  .  (45),  (46),  and  (49)  as 

(50)  v(r,n,  tt)  =  e~^kr  +  •  •  •  > 


(48) 


vg(  r ,  n,  4>)  = 


-  j(kr  +  tt/4) 


V^Trk] 


which  is  one -half  of  the  incident  field  on  the  illuminated  side  of  the 
shadow  boundary. 

The  series  representation  of  vg  given  in  Eq.  (45)  is  valid  every¬ 
where  except  for  the  values 


(51)  4)  =  -fTr  +  2TrnN,  n^2,  N  =  ±  1 ,  +2,  .  .  . 


Near  these  values  the  series  representation  converges  slowly,  and  the 
periodicity  property  of  the  exact  function  v(r,  n,  <p)  can  be  used  to  overcome 
convergence  difficulties  near  these  values.  The  exact  function  v(r,  n,  4>) 
is  periodic  in  2Trn  so  that 
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(52) 


v(r,  n,  4>  +  2TrnN)  =  v(r,n,  <£),  N  =  0,+.l,  +  2,  . 


Therefore,  if  the  series  representation  of  vg  converges  slowly  because 
4>  is  near  one  of  the  values  expressed  by  N  £  0  in  Eq.  (51),  the  periodicity 
property  of  Eq-  (52)  can  be  used  to  represent  v(r,  j))  by  employing  the 
series  representation  near  =  _+  tt(N  =  0).  The  only  case  for  which  all 
boundaries  are  regular  and  the  substitution  of  Eq.  (52)  is  not  necessary 
is  for  the  thin  half-plane,  in  which  case  n  =  2  . 

We  now  examine  the  behavior  of  the  diffraction  coefficient  vg(r,  n,  (j)) 
for  the  case  illustrated  in  Fig.  27  in  which  two  reflected  rays  exist  at 
angles  and  \\)2  •  We  first  examine  the  values  of  4>  =  4*  i^0  for  the  dif¬ 

fracted  rays  corresponding  to  the  directions  of  these  reflected  rays. 


Fig.  27  ,  Case  of  diffraction  by  a  wedge  where 
two  reflected  rays  exist. 
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The  values  of  the  pertinent  angles  are 


(53)  4»i  =  17  -  40  >  4«2  =  (ZniT  - »)  -  i|<0  . 


Thus  the  values  of  <j>  =  +  \\iQ  are 

(54)  i  =  TT  ,  =  (2n7T  -  TT )  . 


For  cpi  =  TT  Eq .  (47)  expresses  the  value  of  the  diffraction  coefficient  and 
thus  the  total  field  at  ijj  =  i^x  is  one-half  the  geometrical  optics  reflected 
field  on  the  illuminated  side  of  the  shadow  boundary,  as  given  by  Eq.  (50). 

However,  for  cj>2  =  2mr  -  tt  and  n  t  2,  Eq.  (47)  is  not  analytic  and 
consequently  the  following  substitution  is  made: 


(55)  v(r,n,cj>)  — ►  v(  r ,  n,  4>  -  2tt  n) 


for  values  of 


(56)  cf>  >  Trn  . 


The  region  over  which  the  subtitution  of  Eq.  (55)  is  made,  as  expressed 
by  Eq .  (56),  is  determined  from  the  symmetry  property  of  v(  r,  n,  cj>);  i .  e  . , 


(57)  v(r,  n,  -q>)  =  v(r,  n,  +4>) 


That  is,  if  the  substitution  of  Eq.  (55)  is  made  for  values  of  4>  <  Trn,  then 
the  substituted  values  of  ^  -  ^Trn  are  closer  to  the  value  of  non-analyticity 
tt  -2mr  than  the  original  values  of  cj)  were  to  the  value  of  non  analytic ity 
2nTT-Tr.  It  is  also  noted  that  the  value  substituted  for  v(r,n,4>2)  is 
v(r,  n,  -tt),  which  corresponds  to  the  correct  value  for  a  shadow  boundary. 

The  equations  given  in  this  Appendix  apply  for  plane -wave  incidence, 
but  they  can  be  used  for  cylindrical-wave  incidence  as  shown  below.  Con¬ 
sider  the  two  situations  shown  in  Fig.  26.  It  is  desired  to  find  the  field  u 

a 
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in  some  direction  £  for  the  wedge  illuminated  by  a  cylindrical  wave  with 
its  source  located  at  (xQJ  a).  By  reciprocity  the  field  ua  is  equal  to  the 
field  u^  which  is  located  at  the  point  ( r  =  xQ,  =  a  )  and  with  a  plane  wave 
incident  from  the  direction  i\iQ  =  £•  The  value  of  u^  is  given  by  Eq.  (47). 
Thus  using  the  property  expressed  in  Eq.  (44)  the  solution  for  diffraction 
of  a  cylindrical  wave  by  a  wedge  becomes 


(58)  ua  =  v(x0>  n,  £+  a)  ±  v(xQ,  n,  a)  . 


The  field  at  infinity  is  given  by  Eq .  (58)  for  a  perfectly  conducting  wedge 
illuminated  by  the  line  source  at  (xQ>  Ol)  . 

If  the  quantity  (akxQ)  is  sufficiently  large  for  Eq.  (48)  to  be  valid, 
then  the  diffraction  pattern  for  cylindrical-wave  incidence  has  the  same 
form  as  that  for  plane-wave  incidence.  In  other  words,  the  diffraction 
for  cylindrical-wave  incidence  is  the  same  as  that  for  plane-wave  inci¬ 
dence  in  regions  sufficiently  removed  from  the  shadow  boundaries  .  The 
region  near  the  shadow  boundary  in  which  Eq.  (47)  must  be  used  increases 
for  decreasing  values  of  xQ  and  may  encompass  all  360  degrees  . 
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APPENDIX  B 

A  LISTING  OF  COMPUTER  PROGRAMS 
USED  IN  THE  PREPARATION  OF  THIS  REPORT 


These  programs  are  coded  in  SCATRAN,  an  automatic  compiler 
language  developed  by  the  Numerical  Computation  Center  of  The  Ohio 
State  University.  SCATRAN  is  similar  to  FORTRAN  IV  in  the  form 
of  arithmetic  expressions  used.  An  experienced  programmer  can  easily 
translate  these  SCATRAN  programs  into  any  desired  language. 
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n  o 


**  * 
**  * 
#*# 

C 


start 


SI  00 

*«* 

00500 

0.1 


RUN 

DUMP  LOWER  CORE 

scatran 

DIFFRACTION  ANALYSIS  OF  bACKSCATTER  bY  A  CIRCULAR  DISC-EDGE 
ON- 

OGlvE  ATTENUATION  AND  DIFFRACTION  COEFFICIENTS  USED- 
wavelength=i • - 

COMPLEX (E 1  *E2*E3«E4«E5*  CEXPL • • DSO • ALPH *  X5 • EC W  *  EC  WT  « ETOT >- 
COMPLEX (E6 >- 
Pl=3. 1415927- 
TP- 2 • #P 1  — 

U 1 =0* 0974- 
U2=0*2 10- 
R 1 3= 1 • /3 • - 
R23=2 • /3 • — 

P 1 6=P 1 /6* - 

PI  1 2-P 1 / 1 2« - 

READ  1 NPUT  «  8 •  ( NK  A )  — 

WR 1 TE OUTPUT •  3  •  <  NK  A )  — 

READ  1 NPuT  «  7 •  ( SFKA*DELKA  )- 

wr i te output •  2  • (sfka*delka  >- 
WR 1 TEOUTPUT  *2- 

E 1 =COS*  (PI  12)— • 1  •  S  1  N  •  <  P  1  12  )  — 

E2- COS • ( P 1 6 ) + • 1 • S 1 N • <  P  1  6  )  - 
E6*C0S«  (Pl4 )  —  •  1 •  S  1  N  • (Pl4  )  — 

DOTHROUGH( SI  00 ) ♦ 1  *  1  *1  • 1 *LE #NKA- 
F1 *1-1- 

FKA*SFKA+F 1#DELKA- 
A -F< A/ TP— 

Fl^FLPFL*( AfRi3)- 

DS0*U1*F1*E1- 

F2= 1  «/FLPFL«  ( A  «R23  >- 

ALPH=U2*F2*E2- 

X3*2  •  *FKA  — 

E  3  =  C0S • (X3)-*1*S1N*(X3)- 

FL  =  P 1  * A- 

FLK*P1*FKA- 

E4  =  C0S#  (FLK)-.  1 .SIN. ( FLK  )  — 

X5*-FL*ALPH  ♦  •U0«- 

E5=CEXPL • ( X5 ) - 

ECW=DSQ*E3*E4*E5- 

ECWT=2**ECW- 

E  SP  «— ( SORT • ( A/2 • )  )/TP- 

ESP*ESP*E6- 

ETOT=ESP-ECWT- 

Y 1  * *REAL*ETOT- 

Y2  =  • 1  MAG  #ETOT - 

FETS0*(Y1*Y1 )+<Y2*Y2>- 

S 1 GMA«4#*FETSQ/<  a*a )- 

WR 1 TEOUTPUT • 2 • < FKA • DSO • ALPH* ETOT • S 1 GMA )- 

WRlTEOUTPUTf 2- 

C0NT1NUE- 

CALLSU0ROUT1NE( >*ENDJO0« ( >- 

ENDPROGRAM ( START ) - 

DATA 


O.  1 
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M*  * 
*  *  * 
#  *# 

c 

c 

c 


start 


SI  00 


*** 

00500 

0. 1 


RUN 

DUMP  LOWER  CORE 

scatran 

DIFFRACTION  ANALYSIS  of  BACKSCATTER  by  a  SPHERE  - 
A  CHECK  ON  REPORT  OSURF  1015-5- 
wavelength= 1 . - 

COMPLEX ( E 1  *E2  »E3*E4 *E5  «CEXPL* • DSQ  «  AlPH  «X5*ECW«ECWT  «ETOT , EPO 
)  - 

PI =3*1415927- 

TP  =  2 • *P 1  - 

U 1 =0.270- 

U2=0.210- 

R1 3= 1 */3.- 

R23=2./3.- 

Pl6=Pl/6.- 

PI 12=PI/12.- 

RADEG=57.29570- 

READ I NPUT  *  8  *  ( NK A ) - 

WR 1 TE OUTPUT • 3 • <NKA>- 

RE AD  I NPUT  *  7 «  ( SFK A « DELK A  )  - 

WR 1 TEOUTPUT  *2«( SFKA « DELK A  )  - 

WR I TEOUTPUT  *2  — 

E 1 =COS. (P112)-.I.S1N.<P112>- 
E  2  =  COS .  (  P  1  6  )  +  •  1  .SlN.  <  P  I  6  )  — 

DOTHROUGH (5100), 1=1 «1 « I.LE.NKA- 
F  1  =  I  —  1  — 

fka=sfka+f i *delka- 
a=fka/tp- 

Fl =FLPFL. ( A tR 1 3  )- 

dso=ui *fi#ei- 

F2=FLPFL. (FKA iRI 3>- 
X 1  - 0 . 750 *P I  - 
F3  =  FLPFL.  (XI  * R23 ) - 
AlPH=F2*F3*E2/( 4. *A )- 
X3=2. *FKA- 

E3=C0S. (  X3 ) - • 1.SIN. ( X3  )  — 

FL=P l*A- 
FLK*PI *FKA- 

E4  =  C0S.<FLK)-.I.SIN.  (FLK  )- 
X5=-FL*ALPH  hI*0«- 
E5=CEXPl • ( X5 ) - 
ECW  =DSQ*E3*E4  *E5- 
ECWT  =2.*ECW- 

EPO  =  -0 • 5* A— •  I  . A/<  4.*FKA  )- 
ETOT=EPO-ECWT- 
Y 1  =  *REAl.ETOT- 
Y2=  •  IMAG.ET0T- 
FETSQ=(Y1*Y1  )  +  <Y2*Y2  )  — 

SIGMA=4.*FETSQ/< a*a>- 

Y3= .REAl.ECWT- 

Y4  =  .  IMAG.ECWT- 

FMECW  =  SQRT.  ( Y3*Y3+Y4*Y4  )- 

PHECW=F AT  AN2 .  <Y4*Y3>- 

PHDEC=RADEG*PHECW- 

WR 1 TEOUTPUT *2 •  ( FKA  t S IGMA • EPO • ET0T . FMECW fPHDEC  >- 
WR 1 TEOUTPUT  *2  — 

C0NT I NUE— 

CALL SUB ROUT  I NE (  )«ENDUOB.(  )- 

ENDPROGRAM ( START ) - 

DATA 


0. 1 
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«*t 


»«* 

c 

c 

c 

c 

start 


F  FMT) 
F  FMT  2 


PUN 

DUMP  LOwtW  COWL 

scatwan 

DIFFRACTION  analysis  of  bACKSCATTER  6Y  A  SPHERE  - 
A  CHECK  ON  REPORT  GSUWF  )8)5-5- 
bISTATIC  SCATTER  OF  A  SPHERL- 
A'AVELENGTH=  l  *  - 

COMPLEX (E 1  «E2 •CEXPL*  *DSQ,  Ai_PH,ECW,ECWT,ETOT,EPO>- 

COMPLEX ( ALPH)  « ALD 1  • AL D2 . E XL • E XL )  *EXL2  )- 

P ! =3* ) 4 1 5927- 

TP=2. «P| - 

U 1 =0.270- 

U2=0*210- 

RI 3=1 •  /  3  •  - 

P23=2./3*- 

P I  6  =P I /6  *  - 

PI )2*PI/>2.- 

RADEG *57 •29578- 

READ  I NPUT • 0  *  ( NK A iNTHETA )  — 

WRITE  OUTPUT , 3* (NKA .NTHETA >- 
READ  1 NPUT  •  7  •  ( SFKA . DELKA . THE TS • DEL TH ) - 
WRITE  OUTPUT  *2 •  < SFKA , DELKA *  THETS*  DELTH )- 
WR 1 TEOUTPUT  *2  — 

E)=C0S*(P1I2)-*1*S1N*(PI12)- 
E2  =  COS • ( P  I  6  )  ♦  • I *SIN* ( P I  6 )  — 

DO  THROUGH ( S I  00 )• Js )  «|  tJ.LE.NTHETA  - 
F  J  =  J— ) - 

THE  TD=THTS+F  J*DEL  TH  — 

THE TAs degrad* the td- 

THET2=P1 -THETA- 

THET  3  =  THE  T2/2  * - 

WRITE  OUTPUT ,FMT) ,< THETD >- 

( 5X , 0*8  I STATIC  ANGLE  «  *.F)5.0>~ 

WRITE  OUTPUT, FMT2- 

( ) 20H  KA  S)GMA  REEPO 

RETOT  1MET0T  FMECW 

)  - 

DOTHROuGH ( SI  00 ) • >  « ) ,) « | *LE.NKA- 
F  I  =  )  -  )  - 

FKA  =  SF K  A  +F  I  *DELKA- 
A=FKA/TP- 
F  K  =  FK  A / A - 
FI =FLPFL« ( A,RI 3)- 
DSO=U) *F ) *E ) — 

F 2  =  FLPFL • ( FK A , R I  3 )  — 

X ) =  0*  750*P I  - 

F 3  =FLPF L •  ( X  I  , R23 )  — 

ALPH  =  F2*F3*E2/ (4* *A  >- 
C0S3-C0S*  <  THET3J- 

EPOr(-)  •  I  *COS. ( THE  T  2 ) / < 2 • *FK A*C0S3*C0$3*C0S3 >  )*0.5*A 

DL  =A  *S I N •  ( THET2 )- 

D2KL»2.*FK*DL- 

ALPH) =ALPH>* I *FK- 

DL  )  -A  *  THE  T  A- 

0L2= ( TP-THETA )*A- 

ALD) = -AlPH I *DL ) - 

ALD2=-ALPHJ#DL2- 

E  XL  =  COS*  (D2KL  )-*  I  .SI  N.  (  D2KL  )- 

EXL) =CEXPL* ( ALD) ) - 

EXL2-CE  XPL  « ( AL02 ) - 

ECWT=-DSO*EXL*<EXLI +EXL2I- 

ETOT=EPO+ECWT- 

YI =*REAL*ETOT- 

Y2* • I  MAG • E T OT - 


IMEPO 

PHDEC 
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non 


FETSQ* < YI #Y l )+<Y2*Y2)- 
S I GMA=4 • #FETSO/ ( A# A ) - 
Y  3= .REAL.ECWT- 
Y4  = • IMAG.ECWT- 
FMECW=SORT. <Y3*Y3+Y4*Y4)- 
PHECW*FATAN2« (Y4*Y3>- 
PHDECsRADEG*PHECW- 

WR 1 TEOUTPUT .2 «  < FK A « S I GM A * EPO * E TOT  * FMECW * PHDEC  >- 
WR I TEOUTPUT  *2- 
CONT I NUE  — 

CALL SU3ROUT I NE <  >*ENQJOB.  (  >- 

ENDPROGRAM ( START ) - 

DATA 


*  JN 

DU*P  LOWER  CORE 
scatran 

BACKSCATTER  FROM  A  PROLATE  SPHEROID- 

SCATTERING  IN  THE  PLANE  90  DEGREES  FROM  THE  MAJOR  AXIS- 
WAVELENGTH*  1  •  - 

COMPLEX ( E I  *E2»E3«E4  *  E5  *  CEXPL • •  DSQ  • AlPH  *X5*ECW*ECWT *ET0T  *EP0 
)- 

START  P  I  a  3  •  I  4  I  5927- 

T P  =  2 • *P I  - 
U 1 =0.270- 
U2  =  0 • 2  I  0- 
R1 3= I • /3 • - 
R23=2*/3.- 
P I 6  =  P I /6*  — 

PI  I 2  =  P I / 1  2. - 
R  ADEG  =  57 • 29578- 
PAT I  0  =  2.- 

READINPuT.Q*  ( NK A  )  - 
WR I TEOUTPUT .3* INK A ) - 
READ  I NPUT . 7.  ( SFkA  • DELKA  )- 
WR I TEOUTPUT*  2*  < SFKA*DELKA  )- 
WR I TEOUTPUT . 2- 

E 1 =COS • (PI  12)-. I.SIN.CPI 12)  — 

E2-COS •  ( P I  6 ) ♦ • I  .SIN.  <  P I  6 )  — 

DOTHROUGH  <  SI  00  >  *  1*1*1*  I.LE.NKA- 
F  I  =  I  -  I  - 

FKA*SFKA4F I#DELKA- 
A -FKA/ TP- 
F I =FLPFL« ( A.R1 3 )- 
DSO  =  U 1 *F 1 #E 1  - 
F2=FLPFL« (FKA*RI 3 )- 
X  I =0*750 *P I  - 
F3  =  FLPFL* ( XI  «  R23 ) - 
ALPH=F  2*F  3*E  2/ <  4  •  *  A  )  - 
X3  =  2*  *FKA- 

E3  =  COS •  <X3)-. I .SIN. <X3>- 
FL  =P 1  * A  — 

FLK*P I *FKA- 

E 4  =  C0S •  ( FLK )- .  I . S IN.  (FLK  )- 

X5=-FL*ALPH  +.1.0.- 

E5=CEXPl. <X5) - 

ECW=DSQ*E3#E4*E5- 

ECWT=2.*ECW- 

F  AM  =  A  #RAT I O- 

EGO*0.5*FAM- 

EPO  =EGO  +  • 1*0.- 

ET  OT  =  -  EPO— EC W  T- 

Y 1 » .REAL .ETOT- 

Y2= • 1  MAG. E  TOT - 

FETSQ= < Y 1 *Y1  )  +  ( Y2*Y2  )  - 

S 1 GMA=4o*FETSQ/ ( A# A )- 

S I GMAtf=2.*TP*FETS0- 

S I GLOG  =  LOG  » (S1GMAW )- 


SI  00 


•  «* 
««  « 
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SI  00 

Y3=.REAl«ECWT- 
Y4=  • lMAG*ECWT- 
FMECWsSORT* (Y3*Y3+Y4*Y4>- 
PHECW  =  FATAN2*  C  Y 4 • Y  3 )  — 

PHDEC=RADEG*PHECW- 

WR I TEOUTPUT. 2* ( FKA . S I GM A ♦ EPO . E TOT . FMECW • PHDEC ) - 
WR I TEOUTPUT  «  2- 

WRI TEOUTPUT *2* ( FK A • S 1 GM A W * S 1 GLOG ) - 
WR1 TEOUTPUT. 2- 
CONT INUE- 

*«* 

00500 

0. 1 

C ALL SUBROUT  I NE (  )= END JOB, (  )- 
ENDPROGRAM ( START ) - 

DATA 

0*1 

*** 

*** 

*** 

RUN 

DUMP  LOWER  CORE 

scatran 

c  backscatter  from  a  PROLATE  SPHEROID  using  geometrical  optics  and 


n  noon 

CREEPING  WAVE  ANALYSIS  - 

FA  =  SEMI  MAJOR  AXIS  - 

FB  a  SEMI  MINOR  AXIS  - 

THT  =  ANGLE  FROM  MAJOR  AXIS  - 

WAVE  =  WAVELENGTH  - 

COMPLEX ( ALPHL .GAMMAL ♦ AECW.PL 1 ♦ PL2 . PHS . DSQ . EC W • ET O T >- 

COMPLEX(El •  E2  • E3 ) - 

COMPLEX (CEXPL* )- 

COMPLEX (CGO.EGOC ) - 

COMPLEX ( EPO ) - 

COMPUTE  LIST  OF  CONSTANTS  - 

START 

PI =  3« 14 15927- 
TP=2.*PI- 
PI 2-PI /2*- 
U I =0*270- 
U2  =  0 • 2 I 0  — 

U2S-0 • 450- 
DELU2=0*050- 
NU2  =  3- 
R 13= I •  /3  •  - 
R23=2 • /3 • — 

R32  =  3 • /2 • - 
Pl4=Pl/4,- 
P134=3**P14— 

P16=PI/6*- 

PI  1 2  *P 1/12*- 
R ADEG=57 • 29578- 
DEGRAD=0*0l 745329- 
El=COS*  <P112)-*I*SIN.(PI12)- 
E2  =  COS*  (  P  1  6  )  +  • 1 *SIN* (  P  l  6  >  — 

E3=COS* (PI34>*,I,S1N,(PI34>- 
FEX=-0 • 5- 
FEXK=  0,5- 

READ  INPUT. 7.  (FA.FB.THTS.DTHT  >- 
WRITE  OUTPUT. 2* (FA.FB.ThTS. DTHT ) — 

READ  INPUT.8.  (NTHT.NWAVE)- 
WR I TEOUTPUT. 3*  (NTHT.NWAVE  >- 
WRITE  OUTPUT. 2- 

fab=fa*fb- 
G3*FLPFL*  <  FAB.R23 )- 

DOTHROUGH ( S 1 00 ) .  I  =  I .  1  .  I  *LE  *NWA VE  - 
READ  INPUT.7. (WAVE)- 
FK  =  TP/WAVE  - 
Gl=FLPFL* ( WAVE ♦ R23 ) — 

G2=FLPFL* ( WAVE.R1 3 )- 
fka=fk*fa- 
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F  FMTI 
F  FMT  5 


F  FMT3 
F  FMT2 

C 


C 


C 


C 


C 


WRITE  OUTPUT.FMT I . (WAVE >  - 
(5X.Q*  WAVELENGTH  =  *.Fl5.8>- 
WR1TE  OUTPUT  • FMT5 •  (FKA  )  — 

( 5X . Q*Fk  A  =  * • F 1 0  #5  )  — 

DO  THROUGH ( S 1 00 ) • KK  =  1  •  I  • KK  *LE • NU2- 
FKK^KK- 

U2  =  U2  S  +  FKK  #  DEL  U2  — 

WRITE  OUTPUT  •  F MT  3  •  (U2 )  — 

(5X.Q*  U2 -  *  .  F  1  0  •  5  ) 

write  output. fmt2- 

(120H  THT  REECW  1MECW  REEGO 

•MEGO  RETOT  I  ME  TOT  SIGMA 

)  - 

COMPUTE  THE  COMPLEX  ATTENUATION  CONSTANT  OF  THE  CREEPING 

WAVE  USING  THE  ELLIPTIC  INTEGRAL  - 
TCW1 =PI 2- 
TCW2=3**P 12- 

CALL  SUBROUT 1NE(PALPH)=F1NCP* (TCW1 *TCW2.0.0l .2iFA.FB.FEX)  - 
ALPHL=U2*E2*PALPH*G3/G2- 

COMPUTE  THE  PATH  LENGTH  USING  THE  ELLIPTIC  INTEGRAL  - 
CALL SUBROUT 1NE (PALPHK ) =F 1 NCP • (TCWl  • TCW2 • 0*0  I  . 2 • FA . FB « FEXK  )- 

fktl=fk*palphk- 

G AMMAL = - ALPHL- • I oFKTL  - 

A  ECW*=  CEXPL  •  (GAMMAL)  - 

DO  THROUGH ( S 1 0 ) • J* 1  .  1  «  J«LE  *NTHT  - 

F  Js  J  - 

THT*THTS^J#OTHT  - 
RTHT=DEGRAD*THT  - 

COMPUTE  POINTS  OF  ATTACHMENT  AND  DEPARTURE  OF  THE  CREEPING 
WAVE  - 

CS1 =COS# (RTHT )  - 

SSI =SIN. (RTHT )  - 
CS2-CS1*CS1  - 
SS2  =  SS 1 *SS I  - 

XP=FA*FA*SSI /SORT • ( F A*F A*SS2+FB*FB*CS2 >- 
XP= • ABS.XP- 
XM=-XP  - 

VSQ  =  F8*FB*(  I  •-<XP*XP/(FA*FA>  )  )  - 
YP=SQRT • ( YSQ )  - 

YM= - YP  - 

ANG 1 =FATAN2*  ( YP.XP  )  - 

TCW 1 =P 1  - ANG I  - 
TCW2=P I +TCW1  - 

COMPUTE  THE  CREEPING  WAVE  .  GEOMETRICAL  OPTICS*  AND  TOTAL 

BACKSCATTERED  FIELDS  - 

FR  =  SORT • ( XP*XP+YP#YP  )- 

SI =FR*  S I N •  ( TCW1 -P I2-RTHT  )- 

S2  =  S 1  - 

FKS1 =FK*S 1 - 

FKS2-FK*S2- 

PL 1  *  I  •  +  • I  *0*- 

PL2* 1 •+• I #0*- 

PHS-PL I +PL2- 

COMPUTE  DSQ  USING  THE  LOCAL  RADIUS  OF  CURVATURE  - 
C  O  SL  =  C  O  S • ( ANG1  >- 
S  I  N’L  =  S  I  N  •  (  ANG1  )  - 
FRLA=FA*SINL- 

frlb=fb*cosl- 

FRLP=FRLA*FRLA+FRLB*FRLB- 
FRL= (FLPFL. (FRLP*R32) )/FAB- 
F 1 = FLPFL. ( FRl . R l 3 ) — 

DSQ=U1*F1*G1*EI- 

ECW=-DSQ*AECW*PHS- 
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F  FMT4 


F  FMT6 

S  1  0 
S  1  00 


S5 


sio 


XGO=FA*FA*CS  1  /SORT.  < F A*F A*CS2+FB*FB*SS2  )  - 
XGO  =  • ABS • XGO  — 

YGO  =  Fd*SORT.  (  1  • —XGO*  XGO / (FA#FA )  >- 
YGO= • ABS* YGO- 
THTG0=FATAN2* (YGO.XGO)- 
FRGO=  SORT •  (XGO*XGO  +  YGO*YGO >- 

rtgo=rtht-thtgo- 
SGO=FPGO*COS* (RTGO ) - 
SSG0=2**FK*SG0- 
CGOrCOS* ( SSGO  1.SIN*  < SSGO )- 
SRGO=FB*FB* SS2+F A #FA*C 52- 
SORGO  =  SOW  T . ( SRGO  )  - 
FACT=0.5*FB*FB*FA/SRGO- 
EPO=-F ACT* (  1  • ♦ • 1  •  1 • / <2*  *FK*SQRGO  )  ) - 

lgoc=e;po*cgo- 

ETOT=EGOC+ECW- 
EMAG1 *  #REAL • E TOT¬ 
EM  AG2*  •  1 MAG*  ET  OT - 
EMSQ=EMAGl*EMAGl*EMAG2*EMAG2- 
S1GMA*2.*TP*EMSQ/<WAVE*WAVE  >- 
WR 1 TEOUTPUT  *  2  «  (THT*ECW*EGOC*ETOT*  SIGMA  >- 
WRITE  OUTPUT  *FMT4 • (SI  )  - 
( 5X  «  Q*  SI  *  *«F10*5>- 
ECW1  =  *REAl*ECW- 
ECW2* • IMAG.ECW- 

FMECW*SQRT • ( EC W 1 *ECW 1 +ECW2*ECW2 ) - 
EGO  1  * .REAL  *EGOC- 
EG02S • 1MAG.EG0C- 

FMEGOaSQRT*  <EG01*EG0 1 +EG02*EG02  >- 
WR 1 TEOUTPUT  «  2  «  ( FMECW • FMEGO )  — 

SIGMAL*LOG* (SIGMA)- 

sigmal«io**s igmal- 

WRITE  0UTPUT«FMT6«  (SIGMAL  )- 
( 5X • Q*S 1 GMAL  *  *«F10*5>- 

write  output  * 2— 

CONTINUE- 
CONT 1NUE- 

CALL  SUBROUT 1NE < >*END JOB* < >- 
SUBROUT 1NE<  SSS )*FINCP.  ( FLL • FUL • ERRR  «NX • FMA • FMB «FEX )  — 
FN*NX- 

DEL* (FUL-FLL )/FN- 
s  SS=0* - 

ERR=0*0 1 *ERRR/FN- 
A*FLL- 

DOTHROUGH ( S40 ) iNNXs  1  •  1  *NNX.LE*NX- 
MXX  =  0  — 

b=a+del- 

ss=o.- 

MX*2- 

DX*0El_/2  •  — 

L  X  *  1  — 

X  *  A  — 

TRANSFERTO( S15)- 
TRAZ=DX*SS- 
MX  *  1  — 

L  X  =  1  - 

dx*del- 

ss=o*- 

LX=LX+1 - 
DX*0*5*DX- 
X  =  A-fDX- 
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S 1 5  DOTHROUGH ( S20 ) *  1X=1 •  1  . IX.LE.MX- 

GS1N=FMA*SIN. <X>- 
GCOS=FMB*COS. <  X  )  - 
GGSS=G  S  l  N*GS  I  N+GCOS»GCOS- 
FSS=FLPFL* (GGSS.FEX  >- 
S 1 8  SS=SS+FSS- 

S20  X=X+2.*DX- 

PROV 1 DED ( LX • E • 1 ) • TRANSFERTO ( S5 > - 
MX=2*MX- 

S  2 1  TRAP  =  0#5*TRAZ  +  0X*SS- 

D  I  F  * • A8S.  ( TRAP-TRAZ )- 
PROV I  DEO (OIF  «GE .D 1 P ) *MXX  =  MXX+1- 
O I P=D I F- 

522  SIMP= (4 •♦TRAP-TRAZ )/3.- 

523  FNCP= ( I 6**S! MP-S1 MZ )/ 1 5.- 

32  4  ER=  *ABS* (  1 .-FNCZ/FNCP )- 

TR  AZ  -  TRAP- 
S  I  MZ-S 1 MP- 


c 

c 

c 

c 


fncz=fncp- 

525  PROV I DED (LX 4 ) • TRANSFERTO < S 1 0  > - 

526  PR0V1DED<MXX.G.4) • TRANSFERTO < S30 > - 

527  PROV I  DEO ( ER • G  *ERR  > « TRANSFER  TO  (SlO>- 

S30  sss«sss+fncp- 

S4  0  A  =  A+DEL- 

S50  CONTlNUE- 

NORMAL  EXIT  - 
END  SU0PROGRAM- 
ORDER  OF  DATA  CARDS  ***- 
CARD  1  F0RMAT*7« < F A • FB « THTS • DTHT )  - 
CARD  2  F0RMAT.8* (NTHT.NWAVE )  - 
NEXT  NWAVE  CAROS  FORM AT t 7 • ( WA VE )  - 
END  PROGRAM (START)  - 
***  DATA 


4*0 


2*0  -2.0 


2  •  0 


0004600003 
4. 18879 
3.59038 
3.14159 
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***  RUN 

***  OUMP  LOWER  CORE 

***  scatran 

c  backscatter  from  a  prolate  SPHEROID  using  geometrical  optics  ano 

CREEPING  WAVE  ANALYSIS  - 

c  E  fielo  parallel  to  plane  of  incidence- 

C  FA  =  SEMI  MAJOR  AXIS  - 

C  F0  »  SEMI  MINOR  AXIS  - 

C  THT  *  ANGLE  FROM  MAJOR  AXIS  - 

C  WAVE  *  WAVELENGTH  - 

COMPLEX ( ALPHL • GAMMAL .  AECW.PL1  • PL2 * PHS • OSQ • EC W . ETO T )- 

COMPLEX  <E I .E2.E3) - 

COMPLEX (CEXPL. >- 

COMPLEX  <  CGO .  EGOC ) - 

COMPLEX(EPO)- 

C  COMPUTE  LIST  OF  CONSTANTS  - 

START  P 1 =3. 1 4 15927- 

TP=2 • *P I - 
PI2=PI/2*- 
U I =0*270- 
R 1 3* I • /3 • — 

R23*2 • /3*  - 
R32=3 • /2 • - 
PI4aPl/4*- 
P I  34  *  3*  *P 14  — 

P I 6=P I /6* - 
PI I 2XP I / I  2. - 
R ADEG*57«  29578- 
0EGRA0»0*0l 745329- 
E I *COS* (PI12)-*I«SIN.(PII2>- 
E2*C0S«  <PI6)+.I.SIN#(PI6>- 
E3*C0S#  <  P  I  34  )  ♦  •  I  •  S  I  N  •  (PI  34  )  — 

FEXK*  0.5- 

READ  INPUT .7. (FA.FB.FKAS.DFKA )- 

WRITE  OUTPUT. 2%  < F A . FB . FK AS . OFK A ) - 

READ  INPUT. 8. ( NK A ) - 

WRITE  OUTPUT. 3. (NKA )- 

REAO  INPUT. 7.  (THT  >- 

WRITE  OUTPUT. 2%  (THT >- 

WRITE  OUTPUT. 2- 

rtht*degrao*tht  - 

F  AB*F  A*FB- 
G3=FLPFL* (FAB.R23)- 

DO  THROUGH(Sl 00).  I  *1  . I .  I  .LE.NKA- 
F  I  =  I- 

FKA»FKAS+(FI-1 . )*DFKA- 
FK  *FK A/F  A- 
WAVE«TP/FK- 
GI*FLPFL* (WAVE.R23)- 
G2=FLPFL* (WAVE.RI3)- 
fka*fk*fa- 

WRITE  OUTPUT . FMT I . ( WAVE )  - 
F  FMT 1  (5X.Q*  WAVELENGTH  •  *.FI5.8>- 

WRITE  OUTPUT. FMT5. (FKA)- 
F  FMT5  (5X.Q#FKA  =  *.Fl0.5>- 


F  FMT  2 

WRITE  OUTPUT. FMT2- 
(120H  THT 

REECW 

IMECW 

REEGO 

imego 

RE  TOT 

IMETOT 

SIGMA 

)- 

65 


c 


COMPUTE  THE  COMPLEX  ATTENUATION  CONSTANT  OF  THE  CREEPING 
WAVE  USING  THE  elliptic  integral  - 

TCW1=P12- 
TCW2=P I  - 

CALL SUBROUT  1 NE (PAlPH ) =F I  NT A •  ( TC W I  *TCW2*5#00*2*FA.FB>- 

alphl=E2*palph/G2- 

ALPHL=2**ALPHL 

C  COMPUTE  THE  PATH  LENGTH  USING  THE  ELLIPTIC  INTEGRAL  - 

C ALL SUBROUT 1NE <PALPHK  )=FINCP.  I TCW 1  • TCW2.  1  .00» 2*FA iFB»FEXK  ) - 

palph<=2.*palphk- 
fktl=fk*p alphk- 
GAMMAL=-ALPHL*»  nfktl  - 
AECW=CEXPL • ( GAMMAl )  - 

C  COMPUTE  POINTS  OF  ATTACHMENT  ANO  DEPARTURE  OF  THE  CREEPING 

WAVE  - 

CS  1  =COS*  ( RTHT  )  - 

SSI =SIN* (RTHT )  - 

CS2=CS 1 *CS 1  - 
SS2*SSl*SSl  - 

XP=FA*FA*SS 1 /SORT • ( F A*F A*SS2 +FB*FB*CS2 >- 
XP= • ABS*  XP- 
XM=-XP  - 

VSQ=FB*FB*( I *-(XP*XP/(FA*FA ) ) )  - 
YP-SORT. < YSQ)  - 
Y  M= -YP  - 

ANG1 =FATAN2*  ( VP *XP  )  - 
TC W I =  P 1 - ANG 1 - 
TCW2-PI ♦  T  C  W  l  - 

C  COMPUTE  THE  CREEPING  WAVE  •  GEOMETRICAL  OPTICS*  ANO  TOTAL 

DACKSCATTERED  F1ELOS  - 
FR  =  SQRT,  ( XP*XP  +  VP*YP  )- 
St  =FR#siNq  (TCW1-P12-RTHT  >- 
S2  =  S 1  - 

fksi *FK*Sl- 
FKS2=FK*S2- 
PL1=1 •+. l.O*- 
PL2  3 1 •  +  • 1  •  0  •  — 

PHS=PL 1 +PL2- 

C  COMPUTE  OSO  USING  THE  LOCAL  RADIUS  OF  CURVATURE  - 

COSL=COS* ( ANG 1 >- 
S 1 NL  =  S I N • ( ANG 1  > - 
FRLA=FA*S1NL- 
FRLB=FB*COSL- 

frlp=frla*frla+frlb*frlb- 

FRL= (FLPFL* ( FRLP  *  R 32 )  )/FAB- 
F 1  -FLPFL  ®  < FRL  *R 13)- 
0SQ=U1*F1*G1*E1- 
ECW=-OSQ*AECW*PHS- 

XGO=FA*FA*CSl /SORT • ( F A *F A *C S2+FB*FB*SS2 ) - 
XGO-*ABS*X  GO  — 

YGO=FB*SQRT* ( 1 • — XGO*XGO/<FA*FA ) )- 
YGO  * • ABS • YGO- 
THT GO  =  F AT  AN2 •  <YGO*XGO)- 
FRGOz  SORT •  ( XGO*XGO+ Y GO*  Y GO )- 
RTGO=RTHT-THTGO- 
SGO  =  FRGO*COS •  ( RT  GO ) - 
SSG0=2*#FK*SG0- 
CGO’COS* ( SSGO ) + • I • S 1 N# ( SSGO )- 
SRG0=FB*FB*SS2+FA*FA*CS2- 
SORGO- SORT •  ( SRGO  )  - 


66 


F  FMT4 


F  FMT6 


SI  O 
SI  00 


S5 


S10 


Si  5 


FACT*0.5*FB*FB*FA/SRG0- 

EPO*-FACT*  C 1 • ♦ • I • 1 • / ( 2 • *FK*SQRGO )  )- 

EGOC*EPO*CGO- 

ETOT  *EGOC+ECW- 

EMAG 1  * • REAL.ETOT- 

EM AG2* •  1 MAG*  ET  OT - 

EMS0=EMAG1*EMAG1+EMAG2*EMAG2- 

S I GMA*2  **TP*EMSO/ ( WAVE* WAVE  > - 

WR I TEOUTPUT • 2  « ( THT  «  EC W  «  EGOC  *  ETOT • S 1 GMA ) - 

write  Output «  fmt4  «  <si )  — 

( 5X • Q*  SI  *  * • F 10.5 )  — 

ECW 1 = • REAL • EC W- 
ECW2=* 1MAG.ECW- 

FmECW*SORT. <ECW1*ECW1 +ECW2*ECW2  >- 
EG01 =*REAL.EGOC- 
EG02  * • 1MAG.EG0C- 

FMEGO*SORT. ( EGO 1 *EGO 1 +EG02*EG02 ) - 
S 1GMAL=L0G* ( S I GMA ) - 
S I GmAL* 1 O*  *S 1 GMAL- 

WRITE  OUTPUT  «FMT6«  <S1GMAL*FMECW*FMEG0>- 

<5X«Q*S1GMAL*  *«Fl0.5«2Xt  Q*FMECW*  * *F 1 5. 7. 2X*0*FMEG0*  *.F15* 
7  )- 

WRITE  OUTPUT «2- 
C  ON  T 1 NUE  — 

CONT 1NUE- 

CALL  SUBROUTINE  <  ) =ENDJOB#  <  >- 

SUBROUT 1NE( SSS)=F1NTA*  ( FLL • FUL • ERRR  «NX«FMA*  FMB ) - 

P  1  -  3 • 1 4 1 5927- 

R23*2*/3*- 

FN=NX- 

DEL* ( FUL-FLL ) /F N- 

sss=o.- 

ERR =0*01 *ERRR/FN- 
A=FLL- 

DOTHROUGH( S4  0 ) «  NNX  *  1  « 1  .NNX.LE.NX- 

MXX  30- 

BsA+OEL- 

SS*0. - 

MX  *2  — 

DX«DEL/2 • — 

LX*  1  - 
X  =  A  — 

TRANSFERTO(SIS)- 
TRAZ*DX*SS- 
MX*  1  - 
LX*  1  — 

DX*DEL- 

ss*o.- 

LX*LX+1- 
DX*0. 5*DX— 

X* A+OX— 

DOTHROUGH ( S20 ) • IX* 1 . 1 • IX.LE.MX- 
S 1 NX  *  S 1 N •  (X)- 
C0SX*C0S. (X)- 
BS 1 N*FMB*S I  NX  — 

BCOS*FM0*COSX- 

asin*fma*sinx- 

ACOS*FMA*COSX- 

AACOS* • ABS. ( ACOS ) - 

FRE*SQRT • ( ACOS*ACOS+BS I N*BS 1  N  )  - 

FRS2* AS  IN* AS  1 N+BCOS*BCOS- 

FRS=SORT. ( FRS2 ) — 

FRP*FRS2*FRS/ < FMA*FMB ) - 
CPS1*-BC0S/FRS2- 
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S  16 
S  1  7 


S10 

S20 


S21 


522 

523 

524 


525 

526 

527 
S  30 
S  40 
550 


S5 


S10 


S IPS  1  * SORT* (  1 *-CPS I *CPSI  >- 

P  S 1 =  F  AT  AN2 •  (SIPS!  .CPS  1  )- 

PROV 1 DED  C PS  1 «L • 0«  0  I  ) .TRANSFER  T0(S16>- 

A0  =  FRE*S1N* (PI -X  >/S IN# (PS  I  >- 

A  0  3 • A0S • (  AO)- 

TRANSFER  TO ( S l 7  )  - 

AO=FMA- 

C  ONT 1 NUE  — 

DXE3AACOS-AO*COS*  (PS!  >- 

BO  =  FMB*SORT.  (  1  • -DXE  *DXE/ ( FMA*FMA )  >- 
FRO=SO*BO/AO— 

RAT  =FRP/FRO- 

XRAT=O#04O*RAT- 

U2*0.20+(  1.4  0/EXPE.  (XRAT)  )- 

FARG=FLPFL* ( FRP.R23 >- 

FSS=U2»FRS/FARG- 

SS  =  SS+F  SS- 

X  =  X  +  2 • *DX  — 

PROV I DED ( LX • E • l ) • TRANSFER TO (S5>- 
MX=2*MX- 

TRAP=0#5*TRAZ+DX*SS- 

D  1  F  =  •  ABS  •  (  TRAP-TRAZ )- 

PROV 1 DED (DIF.GE.D1P) .MXX  =  MXX+ 1  - 

D1P=D1F- 

S 1MP= (4#*TRAP— TRAZ )/3.- 
FNCP* ( 16.*S1MP-S1MZ )/l 5.- 
ER=*ABS* ( 1 •— FNCZ/FNCP )  — 

traz=trap- 

simz=simp- 

FNCZ=FNCP- 

PROV 1 DED (LX *L  *4 ) . TRANSFERTO ( S 1 O  >- 

PROV I DED (MXX#G#4 > . TRANSFERTO ( S30  > - 

PROV I DED ( ER#  G • ERR ) • TRANSFER  TO  (SlO>- 

SSS=SSS+FNCP- 

A  - A-f DEL- 

CONT lNUE- 

NORMAL  EXIT  - 

END  SUBPROGRAM- 

SUBROUTINE  ( SSS  )  3 F 1 NCP • (FLL • FUL . ERRR . NX . FM A . FMB.FEX )- 
FN*  NX  — 

DEL* (FUL-FLL )/FN- 

sss=o 

ERR=0 .0 1 *ERRR/FN- 
A=FLL- 

DOTHROUGH( S40 ) .NNX-1  . 1  . NNX  #LE • NX  — 

MX  X  =0  — 

D  = A+DEL- 

ss=o#- 

MX32- 

DX3  DEL/2. — 

LX  3  1  - 
X  =  A  — 

TRANSFERTO (S 15 >- 

TRAZ=DX*SS- 

MX«  1  - 

LX*  1  - 

DX=DEL- 

ss*o.- 

LX=LX+1 — 

DX=0.5*DX- 
X  * A  +  DX  — 
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S  1 5 

DOTHROUGH  < S20 >  ♦  !  X  =  1  ♦  1  .  1  X • LE .MX - 

GS 1 N*F  MA*  SIN.  <X>- 
GCOS=FM0*COS.  (X  >- 
GGSS=GS ! N*GS 1N+GC0S*GC0S- 
FSS»F|_PFL#  (GGSStFEX  )- 

S10 

S20 

SSsSS+F  SS- 
X*X+2.*0X- 

PROv 1 OEO  <  LX • E • 1  ) ♦  TRANSFERTO  <  S5  > - 

S21 

MX*2*MX- 

TRAP=0.5*TRAZ+DX*SS- 

0  IF*. ABS* ( TRAP-TRAZ )- 

PRO V 1 OED  <0!F.GE*0!P ) tMXX*MXX+i- 

522 

523 

524 

D I P  =  D I F- 

S 1 MP  = (4.* TRAP-TRAZ  ) /3 • - 
FNCP  =  (  16.*5i MP-S IMZJ/15.- 
ER=.ABS. < 1  .-FNC2/FNCP  )- 
TRAZ=TRAP- 
S  I  MZ  =  S 1 MP- 
FNCZ=FNCP- 

525 

526 

527 

53  0 

54  0 
S50 

PROV I OED ( LX • L • 4 )  «  TRANSFERTO < S 1 0 )- 
PRO V 1  DEO ( MXX • G • 4 >  . TRANSFERTO  <S30>- 
PROv ! OEO ( ER . G .ERR ) « TRANSFER  TO  <S10>- 

sss=sss+fncp- 

A= A+OEL- 

CONT I NUE - 

normal  exit  - 

ENO  SUBPROGRAM- 

c  *** 

c 

c 

c 

ORDER  OF  DATA  CARDS  ***- 

OATA  CARD  1  FORM A T , 7 « < F A . FQ . FK AS « DFK A > 
DATA  CARD  2  FORMATiQi (NKA)- 
DATA  CARD  3  FORMAT  «  7  «  ( THT ) - 
END  PROGRAM < START )  - 

**# 

2*0 

00  1  00 
90. 

DATA 

l  .0  0 • 5  0. 1 
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